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Abstract 

We review recent developments in holographic hydrodynamics. We start from very 
basic discussion on hydrodynamic systems and motivate why string theory is an essential 
tool to deal with these systems when they are strongly coupled. The main purpose of 
this review article is to understand different holographic techniques to compute trans- 
port coefficients (first order and higher order) and their corrections in presence of higher 
derivative terms in the bulk Lagrangian. We also mention some open challenges in this 
subject. 
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1 Introduction and motivation 

In this article we shall study the characteristics of hydrodynamic system from string theory 
perspective. 

Hydrodynamics is an effective theory, describing the dynamics of some field theory at large 
distances and time-scales. A fluid system is considered to be a continuous medium. When we 
talk about a small volume element (or "fluid particle") in fluid it has to be remembered that 
the small volume element consists of a large number of molecules (specifically the size of the 
fluid particle is much much greater than the mean free path of the system). The equations 
of hydrodynamics assume that the fluid is in local thermodynamic equilibrium at each point 
in space and time, even though different thermodynamic quantities like fluid velocity (v(x,t)), 
energy (e(x,t)), pressure (p(x,t)), fluid density (p(x,t)) etc may vary. Fluid mechanics applies 
only when the length scales of variation of thermodynamic variables are large compared to 
equilibration length scale of the fluid, namely mean free pathpQ. 

Hydrodynamic description does not follow from the action principle rather it is normally 
formulated in the language of equations of motion. The reason for this is the presence of 
dissipation in thermal media. Due to internal friction called viscosity, a dissipative fluid loses 
its energy over time as it propagates. The fluid without any viscous drag is called ideal fluid. 
In the simplest case, the hydrodynamic equations are just the laws of conservation of energy 
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and momentum (we are considering the fluid without any global charge or current), 

V^T"" = . (1) 

At any space time point the fluid is characterized by d+1 variables (in d dimensions): velocities 
u^(x) (of fluid particles) and its temperature T(x)Q All these variables are functions of space 
and time. Since velocities are time like they follow = — 1, therefore total number of 

variables describing the fluid is d which is equal to the number of equations of motion as in 
equation (lip. 

In hydrodynamics we express T^ v in terms of T(x) and u^(x). Following the standard 
procedure of effective field theories, we expand energy momentum tensor in powers of spatial 
derivatives. As we have already explained local thermodynamic quantities (velocity and tem- 
perature) vary very slowly over space-time therefore their derivatives are very small. Thus it is 
legitimate to express the energy-momentum tensor in powers of derivatives of local quantities 
i.e., ds <C (<9s) 2 , d 2 s <C (<9s) 3 , d 3 s, ■ ■ ■ -C • • • where, s stands for any local quantity. 

At the zero*' 1 order, T^ v is given by the familiar formula for ideal fluids, 

Tp, = (e(x) + p{x))u lx {x)u l/ {x) + p{x)g IMl/ {x) (2) 

where e(x) is energy density and p(x) is pressure^] g^ v is the metric tensor of background 
spacetime. If we also consider the fluid system has conformal invariance (a system which is 
scale invariant) then its stress tensor becomes traceless i.e. = 0. This implies that p = ^rje. 
This is the thermodynamic equation of state. 

At the next order in derivative expansion fluid energy-momentum tensor is given by, 

= (e + p)Ufj,u u + pg^ v - 2o> 



pa p/3 
a v 
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+ (g^V-u, pv = (T + u»u v (3) 



where a^ u is proportional to derivatives of u^(x) and is termed as dissipative part of T^ u (x). a^ v 
is a symmetric tensor. It has been divided into a traceless part and a trace part. The coefficient 
of traceless part is denoted by rj and that of trace part is denoted by (. The coefficient r] is called 
shear viscosity coefficient and £ is called bulk viscosity coefficient. Shear viscosity coefficient 
describes fluid's reaction against applied shear stress where as bulk viscosity coefficient measures 
reaction against volume stress. r\ and ( are called the first order transport coefficients as they 
appear at the first order in derivative expansion of energy-momentum tensor. For conformal 
fluid bulk viscosity coefficient vanishes. Since trace of stress tensor is proportional to ( (from 
equation O)) i.e. , 

T£ ~ CV • u. 

In a similar fashion one can write the expression for stress tensor to the second order in 
derivative expansion. We shall come to the second order hydrodynamics later in section [7} 

Before we go ahead let us first motivate why we use string theory to study different properties 
of fluid dynamics. 



1 Other thermodynamic variables can be expressed as a function of temperature. 

2 Next time onwards we shall drop the functional dependence of local variables i.e. we shall write u v {x) as 
and similarly for other variables. 



3 



1.1 Why string theory ? 

The hydro dynamic behavior of a system is characterized by a set of transport coefficients, like 
shear viscosity, bulk viscosity etc as mentioned above. After Relativistic Heavy Ion Collider 
(RHIC) experiments [2], the study of shear viscosity to entropy density ratio of gauge theory 
plasma has developed lots of attention. The QGP (Quark-Gluon Plasma) produced at RHIC 
behaves like viscous fluid with very small shear viscosity coefficient (near-perfect fluid). Such 
a low ratio of shear viscosity to entropy density is very hard to describe with conventional 
methods. The temperature of the gas of quarks and gluons produced at RHIC is approximately 
170 MeV which is very close to the confinement temperature of QCD. Therefore, at this high 
temperature they are not in the weakly coupled regime of QCD. In fact near the transition 
temperature the gas of quarks and gluons belongs to the non-perturbative realm of QCD. Thus 
usual perturbative gauge theory computations are not applicable to explain RHIC results. On 
the other hand in Lattice gauge theory (a technique to study the properties of strongly coupled 
system) it is difficult to compute real time correlators of energy momentum tensor as the theory 
is formulated on Euclidean lattic^J The AdS/CFT correspondence, at this point, appears to 
be a technically powerful tool to deal with strongly coupled (conformal) field theory in terms 
of weakly coupled (super)-gravity theory in AdS space. Holographic techniques (motivated 
from the AdS / CFT correspondence) to compute hydrodynamic transport coefficients exhibit a 
remarkable quantitative agreement with those arising from numerical fits to RHIC data. 

This motivates us to study properties of hydrodynamic systems from the point of view of 
string theory, in particular, the AdS/CFT correspondence. 

Now, we pause our discussion on fluid system a little and discuss the basic idea and working 
principle of the AdS/CFT correspondence. 

2 The AdS / CFT conjecture and fluid / gravity correspon- 
dence 

In this section we briefly review the AdS/CFT conjecture [3J. The correspondence is itself a 
vast subject and there are lots of good review articles on thi^j Therefore, here we shall briefly 
discuss the important points of this conjecture to introduce ourselves to the notations which 
will be used through out this article. 

The original conjecture states the equivalence between two seemingly unrelated theories: 
Type IIB string theory on AdS$ x S 5 , where both AdS§ and S 5 has radius b, with a five form 
field strength F$, which has integer flux iV over S 5 , and complex string coupling r$ = a + ie~^ 
where a is axion and is dilaton field and M = 4 SYM theory in 4 dimension, with gauge group 
SU(N), Yang-Mills coupling qym and instanton angle 0i (together define a complex coupling 
tym = + p 5 ^ - ) in hs superconformal phase, with gs = a = §^ and 6 4 = ATcgsN(a') 2 . 

3 Plcase look at [3] where the author have found a lattice technique to compute the ratio. But their result is 
far from universality ( rj/s — 
4 See H3 for example 
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The best understood examples of the AdS / CFT correspondence relate the strongly coupled 
dynamics of certain conformal field theories (CFT) to the dynamics of gravitational systems in 
AdS spaces. The description is holographic, as the CFT is living on the boundary of the bulk AdS 
space and thus has one lower dimension than the full bulk spacetime. All the characteristics or 
information of a lower dimensional theory is captured in a theory living in one higher dimension 
and vice versa. The description is also called a duality because when the field theory is strongly 
coupled the gravity theory is weakly coupled and vice versa. 

These boundary theories are quantum field theories of a particular kind. They pertain to 
strongly coupled systems which are difficult if not impossible to study perturbatively. The 
AdS/CFT correspondence seems to be a powerful tool to deal with these strongly coupled 
(conformal) field theory. Since our main target is to understand the properties of strongly 
coupled plasma produced at RHIC0 lets us first understand how such field theory can be dealt 
with via weak gravity computations. 



2.1 The Large 't Hooft coupling (A) Limit 

The large 't Hooft coupling limit corresponds to taking A = gy M N = gsN — > oo, while iV — > oo 
also. This is the weakest form of the conjecture. In this large N limit, the super Yang-Mills 
theory becomes nonperturbative i.e. a strongly coupled super Yang-Mills theory, but the string 
theory side becomes very much tractable. This can be realized as follows: from the relation 
6 4 = 4TrgsN(a') 2 , one can easily check that large A limit corresponds to the radius of curvature 
b being much much larger than the string length, b 2 3> a'. Now, the classical string theory 
effective action is given by, 

C ~ axd'Tl + a 2 a' 2 n {2) + a 3 a' 3 H {3) + ■ ■ (4) 

where, IZ^ is a 2n derivative term constructed out of Ricci tensor or Riemann tensor or Ricci 
scalar. Since the near horizon geometry of D3 brane spacetime is asymptotically AdS, with 
both the AdS 5 and 5*5 has radius of curvature b, the scale of Riemann tensor is set by, 

^111 

7v ~ — — — — rs_/ (5) 

b 2 y/g^Na' y/Xa'' V ; 

Therefore the expansion of effective Lagrangian in powers of a' effectively becomes an expansion 
in powers of (A)~3, 

£~ai(A)^+ 02(A)" 1 + a 3 (A)"i+ • - . (6) 



Thus taking A — > 00 limit corresponds to considering the classical supergravity action instead 
of full classical string theory. 



5 In general, QCD does not have any conformal invariance (ft function is not zero). However, at high 
enough temperature, the QCD plasma is well described by some conformal field theory and thus AdS/CFT 
correspondence can approximate it by a gravity dual. With this assumption we apply this correspondence to 
explain the properties of strongly coupled QGP produced at RHIC. 
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2.2 The Field Operator Mapping 



Another important machinery in the AdS/CFT correspondence is the mapping between bulk 
fields and boundary operators. In [5] Witten proposed that the boundary values of supergravity 
fields acts as a source to the corresponding operator in the field theory side. To be more explicit, 
let us define the string theory partition function as, 



II M ex P [ -,S M 



(7) 



where we collectively denote all string theory fields by Xi- Since AdS is a space with a boundary, 
we need to specify the boundary values for these fields to compute the partition function. Let 
us denote the boundary value of the field Xi by xf ■ Therefore the partition function is, in 
general a function of the boundary values xf' s °f the fields, 



j J [dXi] exp [S[xi\] 



On the other hand in the filed theory side the correlation functions of any operators Oi are 
given by, 

S n Z CFT [{J i }} 



(0^)0^X2) ■■■ O in (x n ) 



SJ il (x 1 )SJ i2 (x 2 ) ■ ■ ■ SJ in (x n ) 



Ji=0 



(9) 



where Z C ft[{Ji}] is given by, 



Z C FT[{Ji}] = Jl[{d<p t ]exp -Si^l + J^jid^lV^F^J^Oiix) 



(10) 



According to [6] there is an one to one correspondence between Xi s and CVs such that 

Z C FT[{J l }] = Z s [xf], with X f = Ji- (11) 



For example, if we want to compute correlation function of boundary stress tensor then the 
corresponding bulk field is metric g^ v . We need to first calculate the bulk partition function 
(string theory/gravity partition function). Then following equation ^ we can find the n-point 
correlation function of boundary operators by taking functional derivative of the bulk partition 
function with respect to the boundary values of the corresponding bulk fields. We shall use 
these results later in this article. 



2.3 Fluid/gravity correspondence 

The power of AdS/CFT is not confined to characterizing only the thermodynamic properties 
of boundary field theories. If we consider a black object with translation invariant horizon, for 
example black D3 brane geometry, one can also discuss hydrodynamics - long wave length devi- 
ation (low frequency fluctuation) from thermal equilibrium. In addition to the thermodynamic 
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quantities the black brane is also characterized by the hydrodynamic parameters like viscos- 
ity, diffusion constant, etc.. The black D3 brane geometry with low energy fluctuations (i.e. 
with hydrodynamic behavior) is dual to some finite temperature gauge theory plasma living on 
boundary with hydrodynamic fluctuations. Therefore studying the hydrodynamic properties of 
strongly coupled gauge theory plasma using the AdS/CFT duality is an interesting subject of 
current research. 

The first attempt to study hydrodynamics via AdS / CFT was in [7] , where authors related 
the shear viscosity coefficient 77 of strongly coupled M = 4 gauge theory plasma in large N 
limit with the absorption cross-section of low energy gravitons by black D3 brane. Other 
hydrodynamic quantities like speed of sound, diffusion coefficients, drag force on quarks etc can 
also be computed in the context of AdS/CFT. 



3 Goal of this article 

In this article we explain how to compute different transport coefficients of strongly coupled 
gauge theory plasma using gauge/gravity correspondence. 

The most famous transport coefficient is shear viscosity coefficient. We plan to study a 
systematic procedure to compute this particular transport coefficient of dual plasma in pres- 
ence of possible stringy correction to the bulk Lagrangian. We first discuss about the famous 
Kubo formula in section |4} In section [5] we show how the shear viscosity coefficient of strongly 
coupled boundary gauge theory plasma depends on horizon value of the effective coupling of 
transverse graviton moving in black brane background. However to completely specify the 
boundary plasma, it is necessary to understand its higher order transport coefficients (coeffi- 
cients which appear in second order in derivative expansion). In section [7] we study the second 
order hydrodynamics. We define a response function in bulk and show that the response func- 
tion flows non-trivially with the radial direction and depends on the full black hole geometry. 
We study the radial flow of this response function (of energy-momentum tensor of dual gauge 
theory) in presence of generic higher derivative terms in bulk Lagrangian and solve these flow 
equations analytically to obtain second order transport coefficients of boundary plasma. We 
consider explicit example to explain this formalism in section [8j Finally we end this article 



with some concluding remarks in section 11 



4 Holographic hydrodynamics 

As we have already mentioned, at present the best available analytical tool one has to study the 
properties of strong coupling hydrodynamics is precisely the AdS/CFT correspondence. In this 
section we briefly discuss how one can apply the knowledge of the correspondence to capture the 
hydrodynamic properties of the plasma. The goal of this section is to introduce Kubo formula 
and use this formula in the context of gauge/gravity duality to compute different transport 
coefficients of strongly coupled plasma. 
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4.1 Kubo formula 



The Kubo formula relates transport coefficients of plasma with their thermal correlators. Let 
us consider the response of the fluid to small metric fluctuation g xy = r\ xy + h xy and consider 
h xy = h xy {t,Xz). The expression for energy momentum tensor in fluid's rest frame is given by 
(from equation ^ and also considering conformal fluid), 

T xy = ph xy + V d, ^ + 0(d 2 ) . (12) 

At the level of linear response theory the one-point function of T xy is linear in h xy and when 
expressed in Fourier space the proportionality constant is simply the thermal retarded correlator 
Gxy,xy (Green's function), 

{Txy) = G xy xy h xy . (13) 



Therefore writing the expression of equation (12) in momentum space we get 



G x % y =p-ioor ] + 0(oo 2 ). (14) 
Thus the shear viscosity coefficient is given bj 



t(«) 
rxy 



7?=lim _^Wy (15) 



Hence, to compute shear viscosity coefficient we need to find the real time retarded Green's 
function of EM tensor. Now, the usual prescription of AdS/CFT to compute boundary corre- 
lator is Euclidean. In principle, some real time Green's function can be obtained by analytic 
continuation of the corresponding Euclidean ones. However, in many cases it is actually very 
difficult to get. In particular the low frequency low momentum limit Green's function (which 
is interesting for hydrodynamics) is difficult to obtain from analytic continuation of Euclidean 
one. The difficulty here is, we need to analytically continue from a discrete set of points in 
Euclidean frequencies (the Matsubra frequencies) u = 2Tcin (n integer) to real values of u. The 
smallest value of the Matsubra frequency is quite large. Hence to get information in small 
co limit that we are interested in, is quite difficult. The authors of [8] have done a detailed 
analysis of this difficulty and have given a prescription to compute the real time correlator. 
This is a well-defined holographic method to compute thermal correlators of boundary theory 
and capture the values of different transport coefficients (for example shear viscosity in first 
order hydrodynamics). However we will follow a different path to compute boundary thermal 
correlator. In the next two sections we try to elaborate on this new technique to compute 
transport coefficients holographically. 

In the next section(s), we use this formula holographically to extract the values of different 
transport coefficients for strongly coupled boundary fluid system. 

6 Similarly other transport coefficients are also captured in G^ R \ We shall consider them in section [7] 



8 



5 First order hydrodynamics 



We plan to discuss the properties of hydrodynamic system order by order in derivative expan- 
sion. As we have already mentioned at the first order we encounter shear viscosity, the only 
non-trivial transport coefficient, for conformal fluid without any other conserved current. This 
section is dedicated to discuss the holographic techniques to compute this transport coefficient 
in generic higher derivative gravity. 

5.1 Hydrodynamic limit in AdS/CFT and Membrane paradigm 

In this subsection we briefly review the proposal described in [9] relating the hydrodynamic 
limit of AdS/CFT to the membrane paradigm. This proposal relates a generic transport coeffi- 
cient of the boundary theory to some geometric quantities evaluated at the black hole horizon 
in the bulk. We will concentrate on the shear viscosity coefficient of the boundary fluid and 
confine ourselves at the level of linear response and low frequency limit of the strongly coupled 
gauge theory. 



Membrane side 

Let us start with classical black hole membrane paradigm, which says that the black hole 
has a fictitious fluid on its horizon. In general, the black hole action can be expressed as, 

^eff = ^out + Ssurf i 

where S out contains integration over space time out side the horizon and S sur f is the boundary 
term on the horizon. Physically, S sur f represents the effect of the horizon fluid on the spacetime. 
Let us consider a general black hole background, 

dS 2 = gMN<ix M dx N = g rr dr 2 + g^dx^dx" (16) 

where M, N runs over d + 1 dimensional bulk spacetime and fi, v runs over d dimensional 
boundary spacetime. This black hole has a horizon at and asymptotic boundary (where the 
dual gauge theory sits) at r b . We assume 5*0(3) invariance in the boundary spatial directions, 
that is all the metric components and the couplings in the theory are only functions of r. We 
consider a small perturbation h xy in the 5*0(3) tensor sector of this metric. We use <j)(r, x^) = h x y 
as the off diagonal component of graviton and in the Fourier space, the perturbation looks like, 

<P(r, k,) = J ^ <f>(r, xn e^ M , k, = (-w, k) . (17) 

The action for this massless perturbation can be written as, 

Sout = - [ d d+1 x^~g -±r (V0) 2 , S surf = [ d d x^ f 1 ^!^ ^ x) (18) 

Jr>r h q{r) Jz V. V^T / 



9 



here, II is the conjugate momentum to for the r— foliation and 7 is the induced metric on 
the horizon. We can interpret S sur f as the effect of the membrane fluid on the spacetime and 
U m b = (^=r) as the "membrane 0— charge". 

Now, following membrane paradigm, the horizon is a regular place for the in-falling observer, 
hence, any physical deformation of the system has to satisfy the in-falling boundary condition. 
The in-falling boundary condition implies that near the horizon rn, 



the deformation should behave as 



~ r 



— rh) iul P, for some constant f3 and 



the solution should be a function of the non singular "Eddington-Finklestein" co-ordinate 
v defined as, 

[Ft 



dv = dt + a — dr. 
V 9tt 

This implies near the horizon r^, the deformation satisfies, 

' Qrr 



(19) 




dtc 



(20) 



The above equation (20) puts constraint on the constant (3 as, 



P 



lg rr (r - r h y 



9tt 



(21) 



With equations ( 17p0 ) and some redefinition of the time, we can also express the membrane 
charge as, 

n m6 = —^-di<P{r h ). (22) 
As per our interpretation, II m b is the response of the membrane fluid induced by and using 



equation (17), in linear response, we define a shear viscosity coefficient for the membrane fluid 

as, 



II mfc = iurj mb (f) 



and thus we get, 



1 



Vmb 



(23) 



(24) 



Boundary Side 

With this much of analysis of the membrane fluid, we concentrate on the boundary side 
where the gauge theory lives. This is a interacting theory at finite temperature and behaves as a 
fluid at sufficiently long length scale or low energy. The real time (Lorentzian signature) finite 
temperature version of AdS/CFT correspondence allows to compute various hydrodynamic 
quantities of this gauge theory at strong coupling by doing some supergravity calculations in 
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the AdS space. Now using the Kubo formula, the shear viscosity of the boundary fluid is given 
in terms of retarded Green's function, the response of graviton to the boundary stress tensor. In 
[8] , the authors have given a simple prescription to compute the boundary correlator using the 
bulk field 0, the off diagonal component of the graviton. Their prescription requires to find a 
solution for the graviton which is in-falling at the horizon and constant at the boundary. Then 
one computes the on-shell action with this solution and the retarded Green's function is related 
to the surface term of the on-shell action at the boundary. Taking 0(fc M ,r) = f(k fM ,r)(f)o(k l j) 
with normalization f{k^r b ) — > 1, 



S 



E 



r=r h ,r b 



d d k 
(2tt)° 



(f>o(k^)G(kfj,,r) 0o(-M 



G R (K 



-99 



lim 2G(ku,r) = lim 2- 
r-*n q[r) 



(25) 



For this profile of the graviton, we can evaluate its conjugate momenta II, it readily gives us 
the relation, 



G R {k, 



lim 



n(A; M ,r) 



r^r b <fi(k^,r) 

Hence, the shear viscosity coefficient rj can be written as, 

n(^,r) 



(26) 



lim lim 



k^or^r b iucpik^^r) 



(27) 



Important point to note is that, in the low frequency limit (k^ — > 0, with II, u<j) fixed), the 
flow of II and u<p in the r— direction are trivial. Hence, we can actually compute the shear 
viscosity coefficient of the boundary fluid at any constant r— slice and its value would be same. 
We compute it at the horizon and get, 



q(r h ) V 9rr9tt 



1 A 

q(r h ) V 



(28) 



Comparing equations (28) and (23), we see that the viscosity coefficient of the boundary fluid is 



related to that of the membrane fluid and more importantly they are given as just the value of 
the inverse effective coupling of the transverse graviton evaluated at the horizon. To emphasize, 



equation (20) plays a crucial role in this equivalence. The AdS/CFT response of the graviton is 



almost same as that of the membrane except that the membrane now has to sit in the boundary. 
The in-falling boundary condition of the graviton field in AdS/CFT is precisely the regularity 
condition given by equation ( 20 ) of the membrane paradigm. In the low frequency limit 



we 



can place a fictitious membrane at each constant r and define the transport coefficient as r](r). 
Since the flow is trivial in this limit, i](r) actually comes out to be a constant, 
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5.2 Example: Two derivative gravity 

In this section explain how to calculate the shear viscosity coefficient of the boundary fluid 
from the effective coupling constant of transverse graviton in Einstein-Hilbert gravity. 

We first fix the background spacetime. We start with the following Einstein-Hilbert action 
in AdS space. 

Here we have taken the radius of the AdS space to be 1. The background spacetime is given 
by the following metric^] 

dv^ 1 1 1 — 

ds 2 = -h t (r)dt 2 + -P— + -dx 2 , Mr) = - -, h r {r) = 4r 2 (l - r 2 ). (30) 

h r (r) r r 

The black hole has horizon at tq = 1 and the temperature is given by, T — - . 
We consider the following metric perturbation, 

9*y = 9 { S + K(r,x) = g<®(l + e$(r,x)) (31) 

where e is an order counting parameter. We consider terms up to order e 2 in the action of 
$(r, x). The action (in momentum space) is given by, (up to some total derivative terms in the 
action^ 

S = / d ^dr (4 O) (r)0'(r, -k)<f>'(r, k) + A { °\r, fc)0(r, k)<j>(r, -*)) (32) 



where, 



<»(r) = ^i, 4\, k) = -^—. (33) 



This can be viewed as an action for minimally coupled scalar field <f>(r, k) with effective coupling 
given by, 

i .4 {0) 

167TG 5 ^/-g(0) g rr 



Therefore according to [5] the effective coupling K e $ calculated at the horizon r gives the shear 
viscosity coefficient of boundary fluid, 

r) = r - l (-2K eS (r )) = — 1 — . (35) 

107rG5 

However, when we add higher derivative terms (effect of string theory) in bulk action the 
situation becomes complicated. In the next two subsections we generalize this idea to higher 
derivative gravity theories. 



7 We are working in a coordinate frame where asymptotic boundary is at r — > 0. 

8 Though throughout this article we have written the four vector k, but in practice we have worked in k — > 
limit. In all the expressions we have dropped the terms proportional to k or its power. 
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5.3 The Effective Action 



Having understood the above procedure of determining the shear viscosity coefficient from the 
effective coupling of transverse graviton it is tempting to generalize this method for any higher 
derivative gravity. The first problem one faces is that the action for transverse graviton no 



more has the canonical form as in equation (32). For generic 'n' derivative gravity theory the 



action can have terms with (and up to) 'n' derivatives of $(r, x). Therefore, from that action 
it is not very clear how to determine the effective coupling. In this section we try to address 
this issue. 



We construct an effective action which is of form given by equation (32) with different 
coefficients capturing higher derivative effects. We determine these two coefficients by claiming 
that the equation of motion for 0(r, k) coming from these two actions (general action and 
effective action) are same up to first order in perturbation expansion (in coefficient of higher 
derivative term). Once the effective action for transverse graviton is obtained in canonical form 
then one can extract the effective coupling from the coefficient of 4>'(r, k)<j)'(r, —k) term in the 
action. Needless to say, our method is perturbatively correct. 

5.3.1 The General Action and Equation of Motion 

Let us start with a generic 'n' derivative term in the action with coefficient /i. We study this 
system perturbatively and all our expressions are valid up to order /i. The action is given by, 

S = — 1 — / d 5 x (R + 12 + fx KM) (36) 
167rG 5 J 

where, IZ^ is any n derivative Lagrangian. ll is the perturbation parameter. The metric in 
general is given by (assuming planar symmetry), 

ds 2 = -(h t (r) + n hf\r))dt 2 + -— + -(1 + ll hf\r))dx 2 (37) 

where h\ n \ hi and are higher derivative corrections to the metric. 

Substituting the background metric with fluctuations in the action given in equation (36) 
(we call it general action or original action) for the scalar field <p(r, k) we get, 

S = y^q- [ tSi^ E 4w( r ' -kH iq \r, k) (38) 

5 J ^ > p,q=0 

where, 0^(r, k) denotes the p th derivative of the field 0(r, k) with respect to r and p + q < n. 
The coefficients A p , q (r, k) in general depends on the coupling constant /i. A p , q with p + q > 3 
are proportional to jx and vanishes in fj, — >■ limit , since the terms (fi^cfi^ with p + q > 3 



appears as an effect of higher derivative terms in equation (36). Up to some total derivative 
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terms, the general action in equation (38) can also be written as, 

1 r d 4 k ^ 
S = 16~G / j2^ dr ^^ r,k ^ iP) ^ r, ~ k ^ iP) ^ r,k ^ n even 

i r d^k (n_1)/2 

~n 77Tv dr E A(r,k)^\r,-k)^\r,k), n odd. 
tG 5 J (27T) 4 f— ' 



16tt( 



p=o 



(39) 



The equation of motion for the scalar field 0(r, fc) is given by, 

where £({0( m )}) is given as 

£({0M}) = ^^ p (r,A;)0^(r,-A;)0W(r,A;) . (41) 
v 

We analyze the general action for the scalar field 0(r, k) and their equation of motion pertur- 
batively and write an effective action for the field <f>(r, k). 

The generic form of the equation of motion (varying the general action) up to order \i is 
given by, 

A (r, k)<p(r, k) - A[(r, k)<p\r, k) - A^r, k)<p"(r, k) = fi £({<P {p) }) + O^ 2 ) (42) 

where J^({0 (p) }) is some linear function of double and higher derivatives of <p(r, k), coming from 
two or higher derivative terms in equation (38). The zero th order (// — > 0) equation of motion 
is given by, 

A ( o\r, k)(f)(r, k) - Al°\r, k)(f)'(r, k) - A?\r, k)(f>"(r, k) = (43) 

where, A^ is the value of A p at jj, — > 0. From this equation we can write 4>"(r, k) in terms of 
4>'(r, k) and (p(r, k) in \\i — >■ limit. 

^f^^fS^' (44) 

Then the full equation of motion can be written in the following way, 

A ( o\r,k)(f>{r,k) -A} 0) (r,k)<f)'(r,k) - Af\r,k)(f>"(r,k) = /i J"(0(r, k),<p'(r, k), <p"(r, k), ...) + C(/i 2 ) .(45) 



Since the right hand side of equation (45) is proportional to /i, we can replace the (j)"(r, k) and 
other higher (greater than 2) derivatives of 0(r, k) by its leading order value in equation (44). 
Therefore up to order /i the equation of motion for is given by, 

Af{r, k)(f)(r, k) - A}°\r, k)<f>'(r, k)— Af\r, k)4>"(r, k) = n J^(0(r, k), 4>'(r, k)) + 0(fi 2 ) 

= n^{r, k) + 7o0(r, k)) + 0(fi 2 ) (46) 

where J-q and T\ are some function of r. This is the perturbative equation of motion for the 
scalar field <p(r, k) obtained from the general action in equation (38). 
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5.3.2 Strategy to Find The Effective Action 

In this subsection we describe the strategy to write an effective action for the field 4>(r, k) whose 



has equation of motion has the form of (32) but with different functions. The prescription is 
following: 

(a) We demand the equation of motion for 0(r, k) obtained from the original action and 
the effective action are same up to order /i. This will fix the coefficients of 2 and 2 terms in 
effective action. 

Let us start with the following form of the effective action. 



1 



dud 3 k 



dr 



(Af\r, k) + fiB^r, k))<f>'(r, -k)4>'(r, k) 



16irG 5 J (2vr) 4 

-(4 0) (r, k) + fiB (r, k))<f>{r, k)<f>(r, -k) 



(47) 



The functions B and B\ are yet to be determined. We determine these functions by claiming 
that the equation of motion for the scalar field <j>(r, k) obtained from this effective action is 



same as equation (46) up to order pi. The equation of motion for (j)(r,k) from the effective 
action is given by, 



4 0) (r, k)cf>(r, k) -X (0) (r, k)<f>'(r, k) - A?(r, k)<j>"(r, k) = p ( fl^r, k) - ^^(r, k) ) <f>'(r, k) 



UAHr^^^-B^k)) <>(r.k) ^, Oi,r) 



(4* 



Therefore comparing with equation (46) we get 



B[(r,k)- 



A} \r, k) 
A?\r, k) 



Bi(r, k) - Ji(r, k) = 



and 



B (r, k) = B 1 (r, k) 



A ( °\r,k) 



-To(r, k) 



The solutions are given by, 



and 



B^r, k) = A?\r, k) \ dr ^ ^ + kA?\t, k) = B x (r, k) + nAf\r, k) 

A\ '(r, k) 



B = B (r, k) + nAf 



(49) 



(50) 



(51) 



(52) 



for some constant k. We need to fix this constant. 

(b) Condition (a) can not fix the overall normalization factor of the effective action. In 
particular we can multiply it by (1 + fiT) (for some constant T) and still get the same equation 
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of motion. Considering this normalization, the effective action is given by, 



1 + n r f dud 3 k 



167rG 5 J (2tt 

|(0) 



-dr 



(A?\r, k) + nBi(r, k))<f/{r, -k)(f)'(r, k) 



+ (<\r, k) + fj,Bo(r, k))<j>(r, k)<j>(r, -k) 



Substituting the values of B's in equations (|51|) and (|52|) we get 
= (I + fi(T + k))S^ + ^ f dr 



Bi(r, k)<f)'(r, -k)(j)'(r, k) + B {r, k)<p{r, -k)<f>(r, k) 



(53) 



(54) 



where is the effective action at fi — > limit. This implies that the integration constant 
k can be absorbed in the overall normalization constant V . Henceforth we will denote this 
combination as V. 

Our prescription is to take V to be zero from the following observation. 

• The shear viscosity coefficient of boundary fluid is related to the imaginary part of re- 
tarded Green function in low frequency limit. The retarded Green function G^ y xy (k) is 
defined in the following way: the on-shell action for graviton can be written as a surface 
term, 

d 4 k 

i>o{k) Gxy,xy(k,r) <f>o(-h) 



S 



(27T) 



(55) 



where (f>o(k) is the boundary value of 0(r, k) and is given by. 



G*{k) = \im2g xy , xy (k,r) 



and shear viscosity coefficient is given 

1 



r->0 
91 



V 



lim 



UJ 



I m Gxy,xy 



(k) 



(computed on — shell) 



(56) 



(57) 



Now it turns out that the imaginary part of this retarded Green function obtained from 
the original action and effective action are same up to the normalization constant F in 
presence of generic higher derivative terms in the bulk action. Therefore it is quite natural 
to take T to be zero as it ensures that starting from the effective action also one can get 
same shear viscosity using Kubo machinery. To show that the above statement is true we 
do not need to know the full solution for 0, in other words to find the difference between 
the two Green functions one does not need to calculate the Green functions explicitly. 
Assuming the following general form of solution for </> 



r 



(5* 



it can be shown generically. For a detailed proof see [TO 



9 To calculate this number one has to know the exact solution, i.e. the form of £ and the value of /3 in equation 
pb. 
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• Because of the canonical form of the effective action, it follows from the argument in [9] 
and the statement above, that the shear viscosity coefficient of boundary fluid is given by 
the horizon value of the effective coupling obtained from the effective action in presence 
of any higher derivative terms in the bulk action. 

(c) After getting the effective action for 0(r, k), the effective coupling is given by, 

K " {r) = IftrfS (59) 

where g rr is the 'rr' component of the inverse perturbed metric and a/— g is the determinant 
of the perturbed metric. Hence the shear viscosity coefficient is given by, 

rj = r^(-2K cB (r = r )) (60) 

where r is the corrected horizon radius. To summaries, we have obtained a well defined 
procedure to find the correction (up to order /i) to the coefficient of shear viscosity of the 
boundary fluid in presence of general higher derivative terms in the action. 

5.4 Membrane fluid in higher derivative gravity 

Let us also consider the effect of the higher derivative terms on the membrane fluid. In higher 



derivative gravity, since the canonical form of the action (32) breaks down, it is not very obvious 



how to define the membrane charge n m b- Instead of the original action if we consider the 



effective action (47) for graviton then it is possible to write the membrane action perturbatively 
and define the membrane charge (II m b) in higher derivative gravity. As if the membrane fluid 
is sensitive to the effective action S eB in higher derivative gravity. 

Following [9] we can write the membrane action and charge in the following way (in mo- 
mentum space) 

«- = / = S^(w^"- fc) ) (61) 

where is the induced metric on the membrane and IT(r, k) is conjugate momentum of 4> 
with respect to r foliation where U(r,k) = {A^ (r, k) + fiBi(r,k))(j)'(r,—k). The membrane 
charge is then given by, 

II mb = ^S- = -k cS (r )y^r dr( p(r, k)\ . (62) 
v— o 



Imposing the in-falling wave boundary condition on 0, it can be shown that the membrane 
charge Il m b is the response of the horizon fluid to the bulk graviton excitation and the membrane 
fluid transport coefficient is given by, 

?7 m b = K eS (r ) . (63) 

Hence, we see that even in higher derivative gravity the shear viscosity coefficient of bound- 
ary fluid is captured by the membrane fluid. 
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Motivated by the fact that like entropy of boundary theory, shear viscosity coefficient is 
also determined in terms of data specified on the black hole horizon, we can think that if we 
know the near-horizon geometry of an asymptotically AdS black hole space time, it would be 
enough to compute the ratio r]/s. In the next section we elaborate this particular issue with 
an example. 

6 Near-horizon analysis of 77/5 

As we know that it is possible to find the entropy of boundary field theory by only specifying 
the complete near horizon geometry of its dual black hole spacetime, one can also ask in the 
same spirit, whether different hydrodynamics coefficients of the boundary plasma can be found 
from the knowledge of bulk near-horizon geometry. There are enough hints from the above 
analysis which seems to point to an affirmative answer to this question. 

In last section we have shown that in the low frequency limit (u — > 0) the shear viscosity 
coefficient rj is related to transport coefficient of membrane fluid, which in turn is given by the 
effective coupling constant of graviton (h xy ) evaluated on the horizon. We have also described 
how to generalize this prescription has been generalized for any higher derivative gravity. All 
these approaches indicate that there can be an "IR" description of transport coefficients rj. 
However, it is not quite clear how other transport coefficients like r n etc. which appear in next 
order in u can be calculated only in terms of horizon data. The derivation of [HI \TU\ was strictly 
valid only in u> — > limit. 

In this section, we observe that only with the knowledge of near horizon geometry of a 
black hole spacetime one can easily calculate the thermodynamic and hydrodynamic quantities, 
namely entropy and shear viscosity of boundary fluid and their ratio. One does not need to 
know the full analytic solutions of Einstein equations. This observation is helpful when bulk 
Lagrangian is very complicated. For example when gravity is coupled to various matter fields 
(gauge field, scalars etc.) in non-trivial way, it may not be always possible to find an analytic 
solution[18j. But if a black hole solution exists then it is possible to find the near horizon 
geometry, i.e. how the metric and other fields behave in near horizon limit. The method is very 
simple. We do not need to solve any differential equation to find the near horizon geometry. 
First we write the field equations. Then take a suitable near horizon ansatz for different fields. 
In general different components of the metric goes like gu ~ o,\{r — r/J + o^r — r/ l ) 2 + - • • , g rr ~ 

^^-(l + 6 2 (^ — rh)-\ )(for non-extremal case), scalar fields behaves as ip ~ iph+(fi(r — r^)H 

and gauge fields goes as A t = q{r — r^) + • • ■ . Substituting these in the field equations we find 
the coefficients consistently order by order in (r — r^). Once we find the near- horizon structure 
of the spacetime we find entropy of the system using Wald's formula, which says we only need 
to calculate the horizon values of some quantities. To calculate shear viscosity we adopt the 
method proposed in last sectionjni [ID] . We write the effective action for transverse graviton 
in black hole near- horizon region and find the coupling constant. For non-extremal black 
hole usually (in presence of higher derivative terms) the effective coupling depends on radial 
coordinate and we need to take r — > limit at the end to find the horizon value of the coupling 
constant. See [11] for detailed discussion and examples. 



18 



However we would like to clarify the following important points : 

Using the AdS/CFT conjecture we are studying the low frequency hydrodynamic properties 
of field theory plasma which has a gravity dual. Our observation is certainly valid for those 
fluids whose dual gravity satisfies the asymptotic AdS boundary condition. The presence of 
a cosmological constant term in the bulk action ensures that the black hole solution must be 
asymptotically AdS. We found near-horizon geometry by solving the background equations of 
motion in presence of the cosmological constant term. This ensures that we are solving the 
near horizon geometry of an asymptotically AdS black hole space time. 

One can also consider Lagrangian with different higher derivative terms. Higher derivative 
corrections typically make energy conditions e.g. weak energy conditions in general relativity 
hard to interpret, as a result of which the underlying solutions could be singular. This means 
that it is not necessarily true that the solutions are asymptotically AdS. However in many 
cases, the higher derivative terms of w- n > kind (contraction of Riemann, Ricci or scalar) can 
be treated perturbatively. Usually if the zeroth order black hole solution is non singular then 
these higher derivative terms (when treated perturbatively) do not create any potential hazard 
in the bulk theory (blowing up of some invariant quantity between horizon and infinity) and 
one gets solutions which are asymptotically AdS [10] , [12] , [13] . For more generic situations, 
if the higher derivative terms are treated exactly, one has to take special care to show that 
there exists a non-singular solution everywhere between the horizon and the boundary which 
is asymptotically AdS. 

The asymptotic AdS boundary condition is true for all the gravity models we studied in this 
chapter. For generic models (where we think our result is most useful) of realistic boundary 
plasmas such as QGP or more generically QCD, one should encounter a dual gravity theory 
with non-trivial matter coupling. A complete analytic solution for such a background is almost 
impossible. For these cases, the existence of the asymptotic AdS solution can be checked 
numerically, and the corresponding near-horizon geometry can be obtained analytically. In 
these scenarios, our observation will be useful to extract information about the coefficient rj for 
the plasma. 

Now we will present a particular model where the above observation will be useful and 
we will compute the shear viscosity of the boundary fluid dual to this system just from the 
knowledge of the near-horizon geometry. 



6.1 The model 

We will focus on a 5-dimensional theory of gravity coupled to a massless scalar and an abelian 
electromagnetic field whose action if 
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S 



EM 



— V / d 5 x^g {R + 12- e a ^F, u F^ - d^ip) + S cs , 

1071-0-5 J 



cs 



16ttG\ 



<fx e^A n F vn F r 



(64) 



(65) 
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We also have to add counterterms to regularize the action — for (64), the counterterms are given in [15] 
(see, also, the nice review |16) and references therein for a more detailed discussion). 
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In this section, we consider a constant moduli potential, V(4>) = 2A = — 12/ Z 2 , and also fix the 
radius of AdS to be I = 1. 

Since the equations of motion for the gauge field simplify, we choose the constant in front 
of the Chern-Simons term to be £/3. The action (64) with various values for a resembles 
(truncated) actions obtained in string compactifications [II] . The coupling ( captures the 
strength of the anomaly of the boundary current. 

The equations of motion for the metric, scalar, and electromagnetic field (F^ 



d^A u 



dyA^) are 



aip(r) 



1 



pa 



d^tpd^tp = 



-g 



-gd*<p) = --ae a ^F, u F^ , 



c 







(66) 
(67) 
(68) 



where we have varied the scalar and the electromagnetic field independently. The Bianchi 
identities for the gauge field are F^x] = 0. 

Since we are interested in a theory for which the Chern-Simons term has a non-trivial 
contribution, we consider the following ansatz for the gauge fieldp] 



A = E(r)dt - —ydx 



B 



xdy — P(r)dz . 



(69) 



The magnetic field, B, is fixed to be constant by the Bianchi identities. Thus, the field strength 
is 



12 



-P'(r) \ 

5 
-B0 




(70) 



/ E'(r) 

-E'{r) 





V p, ( r ) 

where ' denotes derivatives with respect to r. 

Our analysis is on time-independent black hole solutions and so we consider the following 
ansatz for the metric 

_ u(r)dt 2 + e 2V ^ (dx 2 + dy 2 ) + e 2W ^ (dz + Z{r)dtf , (71) 



ds< 



which is compatible with the symmetries of the problem. 

In this case, the horizon is located at (the biggest root of) U{r h ) = and the temperature 
can be easily computed on the Euclidean section — we obtain 



T 



U\r h ) 

47T 



(72) 



11 Plcase note that our P'(r) is the same as P(r) in [19] and our Z{r) in the ansatz of the metric (71) is their 
C{r). 

12 Our convention for the coordinates is (r, f, x, y, z). 
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By using the metric ansatz (71), we can rewrite the Maxwell equations as 

[Q( r yV(r)+W(r)+a V (r)y _ 2(BP'(r) = 0, (73) 

[ e 2V(r)+w(r)+a<p(r) {p ( r ) e ~ 2W ^ P' {r) - Q(r)Z(r))}' - 2(BE'(r) = (74) 

where 

Q{r) = E'{r) + Z{r)P'{r). (75) 



It is easier to work with combinations of Einstein equations rather than using directly (66). 
First, we extract the expressions of second derivatives of the functions that characterize the 
metric in the following way: we obtain W"(r) from (rr)- , V"(r) from (xx)- , and U"(r) from 
(z^)-component of Einstein equations. 

Let us now consider the (tz)-component of Einstein equations in which we replace W"(r), V"(r) 
and U"(r) — we obtain 



2W{r) 



[2V'{r)Z'{r) + 3W'{r)Z'{r) + Z"(r)] - AQ{r)P\r)e aip{r) = 0. (76) 



then considering the 

An important observation, which can be drawn by studying the system of equations ((73 



(76), is that a non-zero magnetic field is not compatible with a constant function Z(r) (and, 
also, P(r)). 

The other (independent) combinations of Einstein equations are obtained as follows: by 
replacing U"(r) in the (rr )-component of Einstein equations we get 

2B 2 e a V (r)-w(r) + 2 Q( r )V^ r) + U{r)[2V'{r) + W \r)f 

+[U(r) (2V'(r) + W'(r))]' + ]- e 2W ^Z\rf - 12 = (77) 

by replacing W"(r) in the (zz)-component of Einstein equations 



-4Q(r)V v(r) + U"(r) + U\r){2V\r) - W\r)\ - 2e 2W{r) Z'i 



r 



,2 



+2U(r) [2V"(r) - 2V'{r)W'{r) + 2V'{rf + <p'(r) 2 } = (78) 
and the last one is, in fact, the (xx)-component of Einstein equations 

e 2W(r)[_ 4B 2 e a V (r) _ s e 4V ^U' (r)V (r) + I2e iv{r) } - 2Q(r) 2 e w{r)+2W{r)+a ^ T) 
+U(r)e 4V{r) [2P'(r) 2 e a ^ r) - 3e w(r) (V"(r) + V'(r)W(r) + 2V\r) 2 )} = 0. (79) 



We also use the ansatz of the metric in the equation of motion for the scalar (67) and so this 
equation becomes 



Due to the non-trivial coupling between the scalar and gauge fields, the equation (67) has a 
non-trivial right hand side. The non-vanishing electromagnetic field may also be understood as 
a source for the scalar field. Thus, the scalar charge is determined by the electric and magnetic 
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charges and so it is not an independent parameter that characterizes the system — this charge 
plays an important role when the asymptotic value of the scalar is not fixed (see [TT]). 

We conclude this section with an observation on the Hamiltonian constraint. A vanishing 
Hamiltonian is a characteristic feature of any theory that is invariant under arbitrary coordinate 
transformations — for our system, we can obtain a first order differential equation by replacing 
W"(r), V"(r), and U"(r) in the (tt)-component of Einstein equations. 

The Hamiltonian constraint, which can be enforced as an initial condition, has the following 
expression: 

2B 2 e a V ( r )-4V(r) + jj^ [_ 2 p/( r )2 e a V (r)-2W(r) + ± V '{r)W'{r) + 2V' {rf - 

V(r) 2 ] + 2Q(r) 2 e^ (r) + 2U'(r)V'(r) + U\r)W\r) + \ e 2W{r) Z\rf - 12 = . (81) 



6.2 Non-extremal Near horizon geometry 

In this section, we will first find the near-horizon geometry of the non-extremal black hole, 
which we will need to compute the shear viscosity to entropy density ratio. In the extremal 
limit, due to the attractor mechanism, the near horizon geometry is universal regardless of 
the asymptotic values of the scalars. We will not present a detailed analysis of the extremal 
branches of solutions here, those can be found in [18]. Here we only present the near-horizon 
geometry of the non-extremal case. 

As in [TJ], we work with a coordinate system in which the solution takes the canonical form 
at the horizon. That is, the field strength Fh and the metric ds 2 H are 



ds 2 H 



q dr A dt + B dx A dy 
r 2 H (dx 2 + dy 2 + dz 2 ), 



p dr A dz, 



(82) 



where q and B are the charge density (of the black brane) and the magnetic field at the horizon, 
respectively. In this way, the gauge freedom is removed and the initial conditions are 



U(r h ) = Z(r h ) = P(r h ) = , V(r h ) = W{r h ) = ln(r h ). 



(83) 



A similar analysis (and numerical solutions) in the presence of the Gauss-Bonnet term but 
without the Chern-Simons term was presented in [20J. 

The generic solutions have a non-degenerate horizon. Near the event horizon, they admit a 



power series expansion of the form (using the definition of the temperature (72 ) in the expression 
off/): 

= ATiT(r - r h ) + u 2 (r - r h ) 2 H , 

= ln(r h ) + ui(r -r h ) + v 2 {r -r h ) 2 H , 

= ln(r h ) + wi(r - r h ) + w 2 {r - r h ) 2 H , 

= zi(r-r ft ) + z 2 (r -r h f H , 



U(r 
V(r 
W(r 
Z(r 
E(r 
P(r 



q(r - r h ) + q x {r - r h f + q 2 {r - r h f + 

p(r - r h ) +pi(r - r h ) 2 H , 

+ Vi(r - r h ) + p 2 (r - r h ) 2 H 



(84) 
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It is important to emphasize that, what is generally called near horizon geometry for a non- 
extremal black hole is just a truncation of the above series expansion. To compute the shear 
viscosity, though, we need also some data at the order (r — r h ) 2 . 

Another observation is that, in principle, one can use a boost transformation in z direction 
to set p = (see [19]). However, the boost transformation is singular at some point outside 
the black hole horizon. In our analysis we keep the value of p non-zero and determine it in 



terms of other horizon data. We will see in section 6^3 that the expressions for the entropy, 
shear viscosity, and their ratio remain unchanged if we set the horizon value of P'(r) to be zero 
(in other words, they do not depend of p). This is expected due to the fact that the physical 
quantities should be invariant under the boost transformations. 



By substituting the ansatz (84) in the field equations, we get the following expressions for 
the coefficients at the order (r — r^) p*| 



Vl 

w x 
P 



2B 2 e aVh + r\ (q 2 e a ^ - 6) 

AB 2 e a Vh _ r 4 ( Aq 2 e a Vh + 3^2 _ 34) 

2AixTr\ : 



&irTr h 
AnTr 4 h 



-2oup h / 

?i = ^r 5 (25 2 [(a 2 + 2)e— +4C 2 ] 

However, in higher derivative gravity theories the only coefficient at the order (r — rh) 2 , 
which we need for viscosity bound computation, is U2- But, for completeness, we present the 
expressions of all the other coefficients that appear at order (r — r^) 2 , in Appendix |pj 

A non-extremal charged scalar black hole is characterized by four independent parameters: 
the mass, electric charge, magnetic field, and also the value of the scalar at the horizon, (fh- In 
this case, the horizon radius (and so the entropy) and the horizon value of the scalar depend of 
the asymptotic boundary data (y?oo)- We will see in the next subsection that this is in contrast 
with the extremal case for which we obtain an attractor behavior at the horizon. 



At first sight, it may seem surprising that the data (85, 246) we need to compute the entropy 
and shear viscosity depend also on zi, a coefficient that we do not compute explicitly. However, 
we will see in the next subsection that the final values of the physical quantities depend in 
fact just on four independent parameters, namely (g, B, r h , (p h ) — we 'trade' the mass for the 
horizon radius and so the independent parameters that completely characterize the black hole 
are the ones mentioned above. 



13 We obtain the results as functions of the coefficients (T, q, B, (fh, rh, and z\) — this will simplify the 
computations of the shear viscosity. 
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6.3 Shear Viscosity to Entropy Density Ratio 

Interestingly, as has already been observed that only with the knowledge of the near horizon 
geometry one can easily calculate the shear viscosity of boundary fluid. One does not need 
to know the full analytic solutions of Einstein equations — this method is especially useful in 
higher derivative AdS gravity theories. 

At first sight, it seems that the dual gravitational mode (4.1) does not generally decouple. 
Interestingly enough, the decoupling occurs when the momentum vanishes and this is what we 
need for the computation of the viscosity in the hydrodynamics limit — a detailed derivation 
of this claim has been provided in Appendix [Ej 



By plugging (31) in the action and keeping the terms at order e (at the first order in e, we 



obtain the equations of motion for gravitons), we get the effective action for the perturbation 



of the form equation (38) 



ii 



Next, we integrate by parts to obtain the bulk action for the graviton in the following form 
(up to some total derivative terms): 

S = [ jL*dr[Ai(r, k)<f>\r, k)<j>'(r, -k) + A (r, k)<f>(r, k)(j)(r, —k)], (86) 

lb7rG 5 J (27T) 4 

where 

i p 2V(r)+W(r),,2 

Mr, k) = -ie 2 ^W)[/(r), A (r, k) = ^ " . (87) 

At this point, it is important to emphasize that there are many total derivatives in this ac- 
tion that do not affect the equations of motion for the graviton. For the computation of the 
imaginary part of the two-point function, the coefficient of the term in the bulk action is 
important. The other total derivatives in the bulk action and the Gibbons-Hawking surface 
term contribution exactly cancel on the boundary [21] . It was also shown [21] (straightforward 
to check for present case) that, in the case of Einstein gravity, the ratio of viscosity and entropy 
density is not affected when the matter fields are minimally coupled. The effective coupling 
is 

1 Axij^k) 1 , , 

167rG 5 y/— gg rr 327rG 5 

and so the viscosity coefficient of the boundary fluid stress tensor is 

V = e 2V ^ +w ^\-2K eff (r h )). (89) 
In this case, the shear viscosity to entropy density ratio turns out to be universal, namely 

V - = h (90) 

S 4:71 

14 The terms that contain the derivatives with respect to the spatial coordinates, x, combine in terms whose 
coefficient is proportional with y 1 . Since we work in the hydrodynamic approximation p = 0, these terms do 
not play any role in our analysis. 
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6.3.1 Four derivative action 



Let us now consider the action (64) supplemented with the most general four-derivative inter- 
actions [23]: 



Shd — Sem + 



a 



16ttG, 



d b x\/—g 



C\R a bcdR a C + C2R a bcdF a F° + C^{F 



+ c A F A + c 5 e abcde A a R bcfg R d /° 



(91) 



Since in supergravity actions the gauge kinetic terms couple to various scalars, it will be 
interesting to understand the role of the moduli in computing the viscosity bound. Unlike [23J , 
our action contains a scalar, if, and the coefficients q depend on the value of (p. This resembles 
the four-derivative supergravity action [24] . 

We treat the higher derivative terms perturbatively and apply the effective action method of 



section 5.3 to compute the shear viscosity coefficient of the boundary fluid. However, to obtain 
the viscosity bound we also need the entropy density. We start by using the Noether charge 
formalism of Wald [25] (see, also, [261 EQ] for a discussion in AdS) to compute the entropy 
density — we will need just the data in Section 3.1 and Appendix B. 

When we add higher derivative corrections to the action, the entropy is no longer given by 
the area law — instead, we use a general formula proposed by Wald 



-2tt 



where e a b is the binormal to the surface H. 



dL 



n 



OR, 



^-ab^-cdi 



abed 



(92) 



By using (92), we obtain the following expression for the entropy density: 



+ 



ot 



r\ [ Cl (3zy h - Au 2 ) - 2c 2 q 2 ] + 0(a' 2 ). 



4G 5 4G 5 

We use the expression of u 2 given in Appendix B to rewrite this expression as 



(93) 



/ 3 

a rr 



4G S G, 




+ 7q 2 e aiph - 6 



c 2 2 

— Q 
2 H 



+ 0{a 



I2\ 



(94) 



As expected, the entropy density depends on four independent parameters, namely (r^, q, B, (ph)- 
Since in Wald formula only the four derivative interactions that involve the curvature tensor 
are important, the entropy only depends on C\ and c 2 (c 5 does not appear because the binormal 
has just rt components and the contribution from this term vanishes). 

To compute the four derivative corrections to the shear viscosity coefficient, we have to find 
the quadratic action for the transverse graviton moving in the background spacetime. As in 
previous section, we consider again the following metric perturbation 



9xy = gi°] + h xy {r,x) 



e$(r, x)], 



(95) 
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where e is an order counting parameter. We find the effective action for transverse graviton 
following the prescription described in previous section. From that effective action one can 
extract the form of effective coupling. 

Evaluating the effective coupling in the near horizon we obtain the shear viscosity coefficient 



7} 
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16ttG r 



a' {c x r\ {8q 2 e aiph + 3z 2 r 2 - 32nTvi - 36nT Wl 



Wu 2 + 48) - B 2 (c 2 - 8c!e a% )) 



(96) 



which can be rewritten as (we use the results in Appendix D ) 
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The ratio of the shear viscosity and entropy density turns out to be 



(97) 
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(98) 



In S — > limit this result matches with [23] 

Let us end up this section with a discussion of the extremal limit. In the absence of the 
moduli, the extremality condition is 2B 2 + r\ (q 2 — 6) = and so the shear viscosity to entropy 
density ratio becomes 



1 
s 



1 a' 

T" + — 
Ait 7r 



B< 



-Ci- 
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+ 0(a 
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(99) 



Therefore, there is a drastic change when the magnetic field is turned on. That is, unlike 
the electrically charged solution studied in [23], the leading correction of rj/s in the extremal 
limit depends on both, c± and c 2 . As expected, in B — > limit our result matches with the one 
of [23]. 

Thus we see, with the knowledge of near-horizon geometry we have been able to compute 
the ratio rj/s of a boundary system which is dual to some non-trivial bulk theory, whose analytic 
solution is not known completely. 

In the next section we will discuss about the transport coefficients, which appear at the 
second order of derivative expansion in energy momentum tensor. Other examples of higher 
derivative correction to first order transport coefficients have been discussed in section |8l 



15 Note that our 'q' is different than l q' of [53]. In [53], q is the physical charge (up to some normalization). 
In our case, the physical charge is 
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7 Second order Hydrodynamics 



In the introduction section we have discussed that hydrodynamics is an effective description 
of field theory in terms of derivative expansions of its local variables. One can write the 
energy-momentum tensor of the system in terms of derivative expansion. We have written 
the form of energy-momentum tensor up to first order in derivative expansion in equation 
(13]). In this section we discuss about the second order expansion of energy- momentum tensor. 
At the second order we encounter many more transport coefficients. Before discussing the 
holographic method to compute these transport coefficients we must mention that second order 
hydrodynamic description is necessary in a relativistic theory, as the first order description is 
not consistent with causality issues [28j EHl [30] . For example, consider a fluid system with some 
conserved current J7 M . The current satisfies the conservation equation and diffusion equation, 

d^J* = 0, Ji= -Dd iP , (100) 

where D is the diffusion constant and p is corresponding charge density. Using these two 
equations one can write, 

d t p - DV 2 p = 0. (101) 

This is a parabolic equation (first order in t and second order in x^'s), and does not satisfy 
causality (see [2H1 [29], [30] for details). In order to restore causality one needs to make the 
diffusion equation hyperbolic. The only possible way to do that is to modify the Fick's law, i.e. 



to introduce a new term in the second equation of (100) 



Ji = -Dd iP - rd t J z . (102) 



The coefficient r is called the relaxation time. Now using equation (102) and conservation 
equation we can write, 

rd 2 p + d t p - DW 2 p = 0. (103) 

This is hyperbolic equation. A propagating solution p = exp[icot + iqz] leads to dispersion 
relation 

- tuj 2 - iu + Dq 2 = 0. (104) 



Therefore in q — > oo limit, the wavefront velocity turns out to be v ~ y D/t, which is consistent 
as long as v < c. Thus, the second order terms in diffusion equation or dispersion relation saves 
the causality problem which appear at the first order hydrodynamics (i.e. if we consider first 
order terms only). 

Like the conserved current J 11 the first order energy momentum tensor also faces the same 
problem. Therefore one has to add higher derivative terms to energy-momentum tensor of a 
relativistic fluid. From the conformal symmetry of the system the form of the second order 
energy momentum tensor can be determined (like first order). In this article we will not discuss 
about the derivation of the second order EM tensor rather, we start with the expression. The 
derivation can be found in 16 [271 E3 EH]. In [27] . the author has developed a Weyl-covariant 
formalism which simplifies the study of conformal hydrodynamics. 



IfiT 



In appendix G.l we have briefed the method of 
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The energy-momentum tensor up to second order in derivative expansion is given by, 



= (e + P)vPu v + Pg^ u - + 6 



fJbV 



0"" = V rn 



+\ 1 a^ x a u ^ x + X 2 a^ x n^ x + A 3 fi (/ \^ A . (105) 
The notations are following, 

A"" = rAt" + , 
D = vPVp (106) 

and for a second rank tensor tensor A^, 

< A »v> = _ A ^A^(A^ + A„ a ) - j-A^A^A^ = A <IXV> . (107) 

It is an easy exercise to check that < A tlv> is transverse and traceless i.e. u^A ll,v> = and 
g^ u A^ v> = 0. Thus T^ u is also transverse and traceless up to second order in derivative 
expansion. 

We see that, we encounter five second order transport coefficients which are named tq, k, 
after For non-conformal fluid there are eight other transport coefficients appear in the 

second order. 

Note that if we consider the hydrodynamic system in flat space, then k term does not 
appear, as the Ricci and Riemann tensor vanishe for flat metric. However, k contributes to the 
two-point Greens function of the stress-energy tensor (please look at equation (108)). The A123 



terms are nonlinear in velocity. They also do not appear when we consider small perturbations 
and thus do not appear in Green's function. The parameter rn has dimension of time and can 
be thought of as the relaxation time. This rn term is enough to get rid of the causality problem. 

The main goal of this section is to discuss a holographic technique to compute the second 
order transport coefficients. If we follow the analysis given in section [4], we find that the 
expression for Green's function up to second order in frequencies is given by, 

G x l xy (u, k) = p - ir/u + r]T n u 2 -%[(d- 3)oo 2 + k 2 } + 0(k 3 ) . (108) 

One important thing to note here (also mentioned earlier) that other three transport coefficients 
do not appear in the expression of Green's function. Therefore we can not compute those 
three transport coefficients calculating retarded Green's function. One needs to adopt different 
holographic methods 



In this section we study the radial evolution of response function defined by equation (113). 
The retarded Green's function of boundary fluid is given by the asymptotic value of this re- 
sponse function. We show that, at low frequencies, the evolution of this response function is 
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independent of the radial direction and hence it can be computed either at the horizon or at the 
boundary. Thus, computing the response function at horizon, the shear- viscosity coefficient, 
which is a first-order transport coefficient of the boundary plasma, can be obtained from the 
characteristics of the membrane fluid. In order to completely specify the boundary plasma, it is 
necessary to understand its higher-order transport coefficients and one needs to move away from 
the low frequency limit. In this case, the response function flows non-trivially with the radial 
direction and depends on the full black hole geometry. Hence, although the boundary plasma 
and the membrane fluid have the same shear-viscosity coefficients, other transport coefficients 
(higher order) can differ and it is not clear how the two are related. One possibility is that the 
higher-order coefficients of the two are related by renormalization-group flow equations. 

We again look at leading Einstein-Hilbert (E-H) action with a negative cosmological constant 
in 4+1 dimensions and study the motion of a transverse graviton in this backgrouncp] The 



action and the solution is given by, equations (29) and (30). 



We study the graviton's fluctuation in this background like equation (31). Here we use 
Fourier transform to work in the momentum space k = {—u, /cj.From this action, we can find 
the conjugate momentum H(r, k^) of the transverse graviton (for r-foliation) and the equation 
of motion from action given in equation (32)[^~ 
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n(r, k„) = 2A { °\r, k)(f)'{r, k), W{r, fc„) - 2 J$\r, k)(f>{r, k) = . (109) 
The on-shell action reduces to the following surface terrnp^ 

s= y^ (0 ( ^ o)(r ' ^ ,(r ' (no) 



=0,1 



Following the AdS/CFT prescription given in [5] (also look at section [5]), the boundary retarded 
Green's function is given as, 

r->0 o (k)0o(-«) 

where, (poik^) is the value of the graviton fluctuation at boundary. Full solution of the graviton 
can be written as 0(r, k^) = (j) (kfj,)F(r, where F(r, k^) goes to identity at the boundary. 
We can rewrite the boundary retarded Green's function as, 



Let us define a response function of the boundary theory a: 

n(r, v 



20 



17 Here we restrict ourselves to five space-time dimensions, but the discussions are quite generic and can be 
extended to arbitrary dimensions. 

18 The factor f/167rGs has been absorbed in a[°^ 

19 Plcase look at [31] for detailed discussion on other boundary terms. 

20 We set the zero frequency part of G to zero, as it gives contact terms 
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where u = ko- This function is defined for all r and k^. Therefor the boundary Green's function 
is given by, 

(114) 



G R {k^) = limiuxik^r) 

i — K) 



We will study the radial evolution of the response function x(k^) from horizon to boundary. 
Differentiating equation (113) and using the equations of motion ( |109 ) we get, 



9tt 



X{k p ,r) 2 T(r) 
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where we define 



E(r) = -2Af\r,k 



9tt 



T(r) = 2A i o\r J k 



9u 

9rr 



(116) 



As mentioned earlier, the flow equation (115) is valid for any value of momentum. This is 



a first order differential equation and we need to specify one boundary condition to solve this 
equation. That naturally comes from the behavior of the equation at the horizon. Demanding 
the solution to be regular at the horizon, we get the following condition, 



E(r)T(r) 



L0 Z 



r=l 



117) 



r=l 



For two derivative gravity this boundary condition implies that 
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E(1)T(1) 
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(118) 



which is independent of k^. Therefore the full momentum response at the horizon corresponds 
to only to the zero momentum limit of boundary response, x(fc M , 1) = xi^n — > 0, r — > 0) p^j 
With this boundary condition, one can integrate out the differential equation (115) from 



horizon to asymptotic boundary and obtain the AdS/CFT response for all momentum k^. In 
particular, it is trivial to see that at (u,ki) — > limit, the flow is trivial 9 r x(/c M ,r) = and 



using the boundary condition in equation (117) we get the first order transport coefficient of 
boundary fluid, i.e. the shear viscosity coefficient coefficients turns out to be r\ = 16 l Gr ■ 

In this section, we will go away from (u, ki — > 0) limit. As we have already mentioned, it is 
possible to integrate the flow equation for any momentum (perturbatively) and we can easily 
find the higher order transport coefficients. The usual Kubo approach to compute these coeffi- 
cients requires the full profile of the transverse graviton in black hole background background 
(solving a second order differential equation), where as, using the flow equation, one can get 
these transport coefficients without explicit knowledge of the graviton's profile. 



21 We choose the negative brunch. The sign of the boundary condition in equation ( 117) depends on the choice 
of coordinate. In our coordinate the boundary is at r 



hence we need to choose the negative branch. 
22 However, in higher derivative gravity we will see that the x(A: M , 1) depends on spatial momentum. 
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7.1 A renormalized response function 



When we solve the flow equation ( 115[ ) to get the boundary response function in general it 



involves divergence at the boundary (r — > 0). These are usual UV divergences and to remove 
them we need to re-normalize the response function properly. 

We follow the holographic renormalization prescription of (32j 05] • As the graviton is mass- 
less, we only need to add the following counterterm to the graviton's action, 



S c = t^tV / d A x^-Me,x)U^{e,x). (119) 
lb7rG 5 J r=s 4 



In momentum space, 



S c = ~F7~7T~ I T^iv 73 ^, k)(g tt u 2 + hk^S, -k). (120) 



647rG 5 J r=s (2tt) 4 
Therefore the renormalized Green's function is given by. 

n(r, fc M ) 



Qr = lim 



o <b(r,kn) 32nG 5 



+ ^r(9 tt " 2 + kk* 



(121) 



However we will study the flow of un-renormalized response function defined in equation 



(113) and we define our renormalized response function as, 

X Ren (r,k„) = x(r,U + - V l — L . (122) 



The counter term will cancel the UV divergences appearing in the expression of x an d we 
will get a finite result at the boundary, i.e. lim^o ^ Ren {r, k^) will be finite. From the above 
analysis, we understand that one can get rid of the UV divergences appearing in the response 
function by following the holographic renormalization technique. But, an important observation 
is, this counter term does not add any finite contribution to the result, it only cancels out the 
divergences. Thus, one can study the flow of the un-renormalized response function and ignore 
the divergences piece to get the finite contribution at the boundary. 



7.2 Second order transport coefficients from flow equation 

In this subsection, we compute the higher order transport coefficients by solving the flow equa- 



tion (115) perturbatively up to order u 2 and kf. This is a non-linear first order differential 
equation. Now, the right hand side of this equation is proportional to u. Hence, to solve x 
to order u 2 , we can replace the leading order solution for x i n the right hand side of equation 
This simplifies the situation a lot as the non-linear equation becomes linear. Now 

Therefore up to order u 2 
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we get, 



+ 0(u 2 ,k 2 



(123) 
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The integration constant for the equation can be fixed form the boundary condition in equation 
(117). Putting the value of the constant, the solution takes the form, 



lim 

r-+o 967rG 5 r 
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3q 2 (r - 1) + u(3u + r(w(log(8) - 3) + 6i)) 
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(124) 



Here we have chosen the four momentum to be k = {u, 0, 0, q}. 

This expression has divergence as r — > (UV divergence) and can be removed by adding 
suitable counter term (as explained in the last section). 

Comparing the finite piece of equation (|124 ) at r — > with the generic expansion of the 



retarded Green's function in equation (108) we get 
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Here T = K These results are in agreement with [281 EHl ED]. Thus, we see that, studying the 
flow equation of the response function we can compute the higher order transport coefficients 
perturbatively. Here, we present the results for the second order transport coefficients, but, in 
general it is possible to go beyond second order. 

At this point, it is not clear why only considering the boundary term from action in equation 



(38) is enough to get the correct results. In usual Kubo approach, one needs to take into 



account the Gibbons-Hawking term also. But as the action in equation (38) has well defined 



variational principle, one does not need to add any Gibbons-Hawking term with it. It has 
been explained in great details in [31] that the boundary terms coming from the original action 
and the corresponding Gibbons-Hawking action are exactly same as the boundary terms coming 
from the action in equation ( 38 ) up to terms proportional to <ft 2 and pure divergence terms. The 
4> 2 terms do not contribute to any transport coefficients _ The divergent terms will get canceled 
by the proper counterterms and hence are not important for finding the transport coefficients. 
Thus it is clear that the effective action will give us the correct transport coefficients for the 
boundary plasma. This observation also holds for higher derivative gravity theory 
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7.3 Higher derivative correction to flow equation 

So far we have discussed the flow equation of two point correlation function of energy-momentum 
tensor of boundary theory whose gravity dual is given by Einstein-Hilbert action (two derivative 
action). But it is not obvious how to generalize this for higher derivative case. The construction 
in the previous section was based on the canonical form of graviton's action given in equation 



(38). AS discussed in section 5.3, in presence of arbitrary higher derivative terms in the bulk, 
the general action for the perturbation h xy does not have the above form as equation (38). 



23 The overall sign depends on the choice of coordinate. 
24 They only contribute to pressure of the boundary theory. 
25 In |31j . it has been proved explicitly for gravity theory. 



32 



Rather it will have more than two derivative (with respect to r) terms like <f)'4>", <fi" 2 etc. In 
presence of these terms it is not possible to bring this action into a canonical form (up to some 
total derivative terms). In this section we consider generic higher derivatives terms in the bulk 
Lagrangian. We follow the prescription of section 5.3.2 to construct an effective action " S e s" 
for transverse graviton in canonical form in presence of generic higher derivative terms in the 
bulk. The effective action and original action give same equation of motion perturbatively in 
the coupling of the higher derivative terms. 

Let us consider a gravity set-up with n derivative action in equation (36). We write an 
effective action for transverse graviton in canonical form, 



S, 
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16ttG, 



d 4 k 



dr 



Af D (r, k)4>'(r, k)(f>'(r, -k) + A% D (r, k)4>(r, k)<f>{r, -k) 



(126) 



with some unknown function A^ D and Aq D . We fix these functions by demanding that the 
equations of motion obtained from the effective action and the original action are same pertur- 
batively in /i (see section 5.3.2). 

The generalized canonical momentum and equation of motion are given by, 

II HD (r, k) = 2Af D {r,k)<j)'{r,k), (n HD (r, k))' = 2A$ B (r, k)<f>{r, k). (127) 

Once we find the effective action for the graviton, we follow the procedure in the previous section 
to obtain the flow equation for the boundary Green's function in generic higher derivative 
gravity. 

The boundary Green's function is given by, 



2Af D (r, k)<p'{r, k)(f)(r,-k) 



lim 

r->0 



which can be written using the definition of canonical momentum as, 

n HD (r, fcj 
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Let us define a response function of the boundary theory in higher derivative theory as, 

n HD (r, k„) 
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Therefore the flow equation is given by, 



where we define 
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(131) 



(132) 
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This is the flow equation for two point correlation function of energy-momentum tensor in 
presence of generic higher derivative term in the bulk action. Therefore integrating this equation 
from horizon to asymptotic boundary one can find the higher derivative correction to the 
transport coefficients at any order in frequency/momentum. 

Like two derivative case here also we need to provide a boundary condition to solve this 
equation. The response function x HD (fc M ,T) should be well-defined at horizon. This implies, 



X HD (^,r) 



£ HD (r)T HD (r) 



r=r h 



(133) 



r=r h 



here the horizon is located at r = r^. 

One important point to mention here is that unlike two derivative gravity where x{^iJ-i r h) 
was independent of k^, X HD (^V») can m general depend on k^. We will see this explicitly in 
the next section. Therefore the full momentum response at the horizon may not be able to 
correspond only to the zero momentum limit of boundary response in higher derivative theory. 

Here we have solved the equation (131) numerically and plotted the function in fig. [I] We 
can see that for non-zero k the horizon value of real part of the response function is different 
than that for zero momentum. 
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Figure 1: Flow of response function for Higher derivative bulk action. 



Like two derivative case, the response function in higher-derivative gravity theory also con- 
tains UV divergences. We need to add proper counter term following the holographic renormal- 
ization procedure to cancel these divergences. A little more thinking also says that in presence 
of any higher-derivative term in the action the structure of the counterterm remains same as 



equation (119). Only the overall normalization constant depends on higher-derivative coupling. 



Thus, similar to the leading gravity, the counterterm in higher derivative gravity also cancels 
out the divergence and does not add any finite contribution to the boundary response function. 
One can study the flow equation of the un-renormalized response function and read off the 
transport coefficients from its finite piece. 

After learning the generic techniques to compute the first and second order transport co- 
efficients of strongly coupled boundary fluid in section [5] and [7| we would like to discuss some 
examples where we can apply these methods to compute them. In the next section we plan to 
discuss some examples of higher derivative gravity, motivated from string theory. 
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8 Examples of 2nd Order Transport Coefficients for Higher 
Derivative Theories 



8.1 Eight derivative correction 

In this section we apply the effective action approach for eight derivative terms in the La- 
grangian. We consider the well known Weyl A term. This term appears in type II string theory. 
Adding this term in the bulk action corresponds to -^72 correction in dual large N theory. The 
five dimensional bulk action is given by, 



S = J d 5 x^g (R+12 + 



where, 



W w = C 



hmnk 



a 



1'innq^h ° rsk <~ g 



sk 



and the weyl tensors C abcd are given by, 

C abcd = R abcd + \{g ad R cb + g bc R ad - g ac R db - g bd R c 
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{ga C gbd — g a dgcb)R ■ 



The background metric is given by [33] (with horizon at Tq — 1), 



4(1 - r 2 )r 2 



(1 + 45/ir 6 - 75/ir 4 - 75/ir 2 ) dt 2 



(l - 285/xr 6 + 75/ir 4 + 75/jr 2 ) dr 2 + -dx 2 



The temperature of this black hole is given by T = - (1 + 15//) . 
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8.2 The General Action 



Putting the perturbed metric in equation (134) we get the general action for the scalar field 



r, k). As mentioned in section 5.3, we see that in presence of higher derivative terms the 



general action does not have canonical form. The action for is given by 

1 r d 4 k 
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Expressions for Af can be found in [10] . Up to some total derivative terms this action can be 
written as, 
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where, 
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(139) 



With this form of action it is not possible to define the effective coupling constant of transverse 
graviton or a response function (following equations (34) and (113) respectively). Therefore we 
find the effective action for transverse graviton <p in presence of W 4 term in the action. 



8.3 The Effective Action and Shear Viscosity 



Following the general discussion in section 5.3.2 we find the effective action for graviton in 
presence of W 1 term. We write the effective action for the scalar field in the following way, 
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;i + r» r d 4 k 
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(140) 



The functions £>g and /S^ are given by. 



«?>.*) = -" 2(663 ^ r 5 2 7 ! r |) +75r2> . TM) = (r2 - 1)(129r6 r +141r4 - 75r2> ■ (141) 

We set the normalization constant T = (see [TU] for detailed discussion). Then the effective 



coupling constant is given by equation (59), 

1 Af\r,k) +/i BY (r,k) 1 
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Therefore the shear viscosity is given by, 

77 = r " f (-2ir cff (r )) 

and shear viscosity to entropy density ratio 

77 1 
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where entropy density s is given by s = ^ (1 + 60 /j) [33]. These results agree with the one 
in the literature. 
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8.4 String theory correction to flow equation 

Now let us compute string theory corrections to the second order transport coefficients following 



the general discussion given in section 7.3 



From the effective action for transverse graviton, computed in equation (140), the flow 
equation is given by, 
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where, 
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where, 
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The explicit expressions for S and T can be obtained by using A^ 4 and A^ 4 . From the 
regularity of x W {ku, r ) at horizon we get, 
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Here, we see that unlike the two- derivative gravity, the horizon value of the response function 
depends on spatial momenta q. With this boundary condition we solve the flow equation up to 
order u 2 and q 2 (ignoring 0(uq 2 ) term). Here we write the final resultp 
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Comparing this result with equation ( 108 ) we get 
V 



1 + 135// + C(/i 2 ), k = \ (1 - 145//) + C(/t 2 ) 
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These results are in agreement with [31] who applied usual Kubo formula to obtain these results. 
The agreement provides a non-trivial check to this approach of obtaining higher order transport 
coefficients from the flow equation (115). 



2(i 



k = {uj, 0, 0, q} and we ignore the UV divergence piece. 
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8.5 Four derivative correction 



Next, we will concentrate on the generic four derivative corrections to Einstein-Hilbert action. 
These terms arise in the effective action for the heterotic string theory. In fact, the complete 
super-symmetrized R 2 correction to effective Heterotic string theory is known and one way to 
obtain it is the super-symmetrization of the Lorentz Chern-Simons terms [35j EE] . This terms 
also arises in the context of Type IIB string theory [T3l 137] . where the theory is on AdS$ x X 5 , 
the compact space X 5 being S 5 /Z 2 . The dual theory is Af = 2Sp(N) gauge theory with 4 
fundamental and 1 antisymmetric traceless hyper-multiplets. This super-conformal theories 
arises in the context of iVD3-branes sitting inside 8 D7-branes coincident on an orientifold 
7-plane. In this case, generic four derivative R 2 correction comes form the DBI action of the 
branes. 

Here we compute the generic four derivative correction to the second order transport coef- 
ficients, the relaxation time r n and k. We can choose the coefficients of the higher derivative 
terms to be the four-dimensional Euler density and get pure Gauss-Bonnet correction to these 
coefficients. 

The action 



R + 12 + a'[ fcR 2 + (3 2 R, vpa R^ + fcR^IT 



(151) 



In particular for Gauss-Bonnet correction, (3i = 1,02= l,/?3 = — 4. One can get rid of the 
Ricci 2 and Scalar 2 terms by a field redefinition and therefore all physical quantities should 
depend on the coefficient 0% only. Here, we prefer to work with the generic case as it would be 
easier for us the get the results for pure Gauss-Bonnet combination at every step. 

The background solution is given by [38J, 



ds 2 = f{r)dt 2 + ^dr 2 + ^-dx 2 (152) 



where /(r) and g(r) are given by, 



f(r) = r - l --2r{r 2 (153) 



and 



r 



y{ ' l-r 2 3(r 2 -l) K ' 

This is the background metric corrected up to order a'. We have fixed the integration constant 
such that the boundary metric is Minkowskian and the horizon is located at r = 1. The 
temperature of the black brane is given by, 

T = 1 + 10ft - 5ft + 2gy. (155) 

7T 37T 
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Similar to the Weyl A case, we can write the following effective action for this model, 
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where, Af B and Aq B are given in appendix [cj Now, it is straightforward to write the corre- 
sponding flow equation (115) in this case, 
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where we define 
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Now, the boundary condition (117) takes the following form, 
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As mentioned earlier, we see that even in this case, the boundary condition depends on spatial 
momenta q through the coefficient {3%. With this boundary condition, one can solve the flow 
equation (157) and the solution is given by, 
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From this expression we get the following transport coefficients, 
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This matches with results in [T3| |39| HQ]. Now, we find the higher order coefficients, 
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As we can see, the physical quantities rj/s, k, t^T only depend on the coefficient 02- In 
particular to Gauss-Bonnet combination, the corrections are, 
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1 - 10a') + 0(a 
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2-ln2 11 , m 



(164) 
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Thus, we show that studying the flow of the response function (constructed from the effective 
action), we can find out all higher order transport coefficients systematically. 



8.5.1 Exact result for Gauss-Bonnet black hole 

As we have done the above computation perturbatively, the above expressions are valid only at 
order a'. But, one can consider the Gauss-Bonnet term exactly in coupling. For pure Gauss- 
Bonnet combination the equations of motion remain second order differential equation and 
hence it is easy to solve exactly to find the background space-time. We solve the flow equation 
in this background exactly in coupling constant, and find the exact expressions for relaxation 
time r n and k. In this section we briefly outline the result. 
The action and the solution is given by, 
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In this coordinate the boundary metric is rj. We also consider only the '— ' branch of f± 
which corresponds to a non-singular black hole solution with non-degenerate horizon. 
The black hole temperature is given by, 



T 



r 



(169) 



With exact GB, the effective action for fluctuation has a canonical form. Therefore we 
derive the flow equation for the response function (as we did in section [7]) and solving this 
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equation we get, 

1 /-. iX X _ 2Agb (8 Xgb — 1) 

77 " 1GkG 5 {1 AAsbh K ~ (1 - v/T^V) (4A 9fe - 1) 
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(170) 



One can easily check that up to first order in X gb , the results in (170) reduces to the one in 



(164). In jH] the authors obtained the relation between second order transport coefficients and 



Xgb numerically, however we are able to present the result exactly. 

Violation of bound on shear viscosity to entropy density ratio 

In our units the entropy density s turns out to be s = j^- and hence shear viscosity to 
entropy density ratio turns out to be 

2 = ^(1-4^), (171) 

which violates the famous KSS bound for X gb > jlHl [13]. In (12] the authors conjectured 
a bound on shear viscosity to entropy density ratio for gauge theory plasmas which have a 
holographic dual, i.e., 

- > ^- (172) 

S 47T 

In [10] the authors argue that when X gb > ^ the theory violates microcausality and is in- 
consistent. Therefore, for (3+l)-dimensional CFT duals of (4+l)-dimensional Gauss-Bonnet 
gravity, consistency of the theory requires 

V ->±--. (173) 
s ~ 4vr 25 v ; 

As we have mentioned in introduction that the flow equation is a first order non-linear 
differential equation but one can reduce this equation to a second order linear differential 
equation. This second order differential equation is related to the equation of motion for 
transverse graviton (gauge invariant excitations). Therefore we can use this equation to study 
causality violation in Gauss-Bonnet gravity. In [4"T| |4"3] it was found that to preserve causality 
of a conformal fluid there exists a bound on second order transport coefficients, 

r^T - 2^ > . (174) 

s 

In Fig|2]we plot for our result and find the following bound on X g b which is in agreement 

with [41] . 

- 0.711 < X gb < 0.113 . (175) 



41 



2<1 

Tr, 



O.Od 



0.04 



0.02 



0.00 





10.113 







-O.S -0.6 -0.4 -0.2 0.0 0.2 



Figure 2: Bound on Ag&. 

9 Flow equation for charged black holes 

Electrically charged black holes in five dimensions have drawn a lot of interest in the context 
of AdS/CFT. The electric charge of these black holes are mapped to the global R-charge of the 
dual field theory. Because of the presence of the electric charges, the thermodynamics and the 
phase structure of these black holes are rather complicated and also interesting at the same 
time. There have been a lot of study of thermodynamics and phase transitions of these charged 
black hole with different horizon topologies (see jH] and references therein). 

The goal of the present section is to apply the AdS/CFT correspondence to understand how 
non-vanishing chemical potentials effect the hydrodynamic behavior of strongly coupled gauge 
theories. We study the second order hydrodynamics in two cases: (a) Generic R-charge black 
holes and (b) Charge black holes in higher-derivative gravity. 

9.1 R-charged black holes 

We consider a conformal field theory with conserved charge (density) in addition to energy and 
momentum. This is especially an interesting extension of the hydrodynamics of the uncharged 
fluids. 

The second order hydrodynamics of charged fluid has been studied in (151 ;46j . They consider 
Reissner-Nordstrom black hole in five dimensions and found the effect of chemical potential on 
second order transport coefficients in some limits of chemical potential. One important outcome 
of their analysis was to find a new non-dissipative contribution to the charge current. How- 
ever, we consider generic parity preserving .R-charged black holes with three (unequal) charges 
(chemical potentials) and find the exact expressions for second order transport coefficients in 
presence of three chemical potentials. As we have mentioned in the introduction that solving 
the flow equation (of retarded Green's function of energy momentum tensor) we can only find 
two second order transport coefficients whereas in [45] H6] all other second order transport 
coefficients have been reported. 

We consider R-charged black holes in five dimensions. A consistent truncation of M = 8, 
D = 5 gauged supergravity with SO (6) Yang-Mills gauge group, which can be obtained by 
S 5 reduction of type I IB supergravity, gives rise to H — 2, D — 5 gauge supergravity with 
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U(l) 3 gauge group. The same theory can also be obtained by compactifying eleven dimensional 
supergravity, low energy theory of M theory, on a Calabi-Yau three folds. The bosonic part of 
the action of M = 2, D = 5 gauged supergravity is given by [4"4"j . We follow the notation of 



^sugra t r j d X y/ <7 

107T(jr5 



-sugra 
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R + V(X) - l -G u (X)Fl u F^ J - GuWd.Xjd^Xj 



+ 77rV / d * x ^ vp<T1 ^F up F al (176) 
lo7rG 5 J 

where, X z, s are three real scalar fields, subject to the constraint X 1 X 2 X 3 = 1. F 7 's, which are 
field strengths of three Abelian gauge fields (I,J=1,2,3), and the scalar potential V(X) is given 
by, 

F% = 2d { »A% Gu = ^diag [(X 1 )- 2 , (X 2 )" 2 , (X 3 )' 2 ] , V(X) = 2^i- (177) 

The three-charge non-extremal STU solution is specified by the following background values 
of the metric 

fk = k-^+r 2 n, Hi = l + ^, H = H 1 H 2 H 3 , (179) 
as well as the scalar and the gauge fields 



■"i y Hi 

The parameter k determines the spatial curvature of dVl\ k : k — 1 corresponds to the metric on 
the three-sphere of unit radius, k — - to the metric on R 3 . as the Hydrodynamic approximation 
is valid only in the case of a translational-invariant horizon, in our case we set k = and 

dfi^o -»■ (dx 2 + dy 2 + dz 2 ) . 

Replacing the radial coordinate r — > ro/\/r, where ro is the largest root of the equation f(r) = 0, 
the background solution in this new coordinate is given by, 

dsl = -n- 2 / 3 ^^fdt 2 + n 1 / 3 ^dr 2 + H 1 / 3 ^^(dx 2 + dy 2 + dz 2 ) (181) 



f(r) = -H(r)-r 2 f](l + «:,), #, = 1 + mr , «< = § ' (182) 
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where k'^s are chemical potentials and 

The scalar fields and the gauge fields are given by 



X 1 



Hi(u) 
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where, 
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The Hawking temperature of the background (181 ) is given by 

T _ 2 + Kj + K 2 + K 3 - KiK 2 K 3 
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We perturb the xy component of background metric and the action for transverse graviton 
is given by, 
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Therefore the flow equation is given by, 
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Solving this equation perturbatively in u and q we get 
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0(qcu 2 ,ujq 2 ,q 3 ,uj 3 ) (191) 



where, 
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(192) 



Computing the response function at the boundary (throwing away the divergent piece) we 
get the following transport coefficients, 
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and 
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These are the new results in this paper. It is easy to check that for Ki — > limit we recover the 
results in section [3 

To complete the discussion on the second order transport coefficients for i?-charged black 
holes one should find the flow of Green's functions for two point correlation functions of R- 



currents. As we will mention in section 10 that in presence of finite charges (or chemical 
potentials) it is very hard to solve the Riccati equation even perturbatively in u and q. We find 
it very difficult to get any analytic solution for i?-current Green's function. However we consider 



a simple model in section 10 and study the flow of i?-current Green's function numerically. 



9.2 Charged black holes in higher derivative gravity 

In this section, we will study five-dimensional gravity in presence of a negative cosmological 
constant and coupled to U(l) gauge field. The model has been studied in [23" ] H8f l4"9~] . the action 
is given as, 
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Here, F 2 



F ab F ab ,F* 



F ab F F cd F , and the AdS radius is set to unity. The action in- 
cludes the Chern-Simon term and also a generic set of four derivative terms. All the four 
derivative terms will be treated perturbatively in our computation and here a' << 1 is the 
perturbation parameter. In [23], it was shown that, within perturbative approach, after using 
field-redefinition, this is the most generic four derivative action that one can write down. In 
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this section, we will closely follow their work. The background metric and the gauge field in 
presence of these higher derivative terms have the following form, 



-r 2 f(r)dt 2 + 



r 2 g(r) 



dr 2 + r 2 (dx 2 + dy 2 + dz 2 ), A = h(r)dt, 



(196) 



where, 



f(r) = f (r)(l + a'F(r)), g(r) = / (r)(l + a\F(r) + G(r))), 
h(r) = h (r) + a'H(r). 



(197) 



Here fo{ r ),9o{ r ) and ho(r) are the solution of the background in absence of the higher-derivative 
terms in the action and they are given as, 
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Here, Q is related to the physical charge of the system and r is the position of the horizon. 
From (197), it is clear that even in presence of the higher-derivative terms, the horizon remains 
at To- The higher-derivative corrections to this background are given by the functions F(r), G(r) 
and H(r). The form of these functions are given in [23]. We would not write those expressions 
and refer the reader to that paper. 

Using the flow equation, we will study the higher order transport coefficient of the plasma 
theory dual to this gravity model. For this, we will write the effective action for the metric 
fluctuation in (31), as we have done in previous sections, 
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where, Af B and Aq B are given in appendix [fJ The corresponding the flow equation (115) for 
this case with the coefficients A± B and Aq B is, 
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We solve this flow equation to find the effect of higher derivative terms and chemical poten- 
tial (or charge) on transport coefficient. Here, we present our result for small Q only, though 
it is possible to find the results for any Q. 



The boundary condition (117) will take the following form 27 
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Note that we are working in a different coordinate where r — > oo is the boundary, therefore we choose the 
positive branch of the boundary condition (117) 



46 



In this case, the horizon value of the response function is independent of momenta. With this 
boundary condition, we can solve the flow equation (200), and the solution is given by, 
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It is easy to read off the transport coefficients from this expressions. 
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where, the temperature T of the system is given by, 
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We see that the first order as well as the second order transport coefficients coming from 
retarded Green's function of energy momentum tensoip^] only depends on two coefficients ci, c 2 . 
This feature was observed in [23] for entropy density s and first-order transport coefficients rj. 
The coefficients C3, C4, which parameterize couplings in the four point function of the dual U(l) 
current does not play any role in these hydrodynamic coefficients. They should be important 
for the computation of conductivity, which comes from the Green's function of the boundary 
i?-current. Two other coefficients (, c§ also do not appear in the expressions. One can find a 
magnetic brane solution of the action (195) like [19J. In that case it would be interesting to 



find the effect of magnetic field on transport coefficients. 

So far we have discussed the flow equation for energy-momentum tensor. However, in 
general the fluid can have other conserved currents like R current, if the theory has a global 
R symmetry. Therefore, in the same spirit, one can also discuss the flow of R current in the 
context of fluid/gravity correspondence. This is the final topic of our discussion. 



28 It would be interesting to study the flow of retarded Green's function for boundary i?-current. We found 
it to be difficult to get any analytical solution for response function in presence of finite chemical potential and 
higher derivative terms. However, it would be nice to know the higher derivative corrections to other second 
order transport coefficients appear in R current [35] I46j . 
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10 Flow of retarded Green's function of boundary R cur- 
rent 



Finally, in this section, we study the flow of retarded Green's function of boundary i?-current, 

Gfm =-i f dtd 3 xe ik - x ([Ji{x), J^O)]) (206) 



where J M (x) is the CFT current dual to a bulk gauge field A^. 

In hydrodynamic approximation one can express the current in powers of boundary deriva- 
tives. Up to first order in derivative expansion it has the following form, 

■h = -%P?d a ^ + Ql v + 0{d 2 ) (207) 

where, k and Q are two first order transport coefficients, fi is chemical potential, T is temper- 
ature and 

l» = efi Ua dp Uj . (208) 

The expression of J M up to second order in derivative expansion can be found in |l5j H6] . From 
conformal invariance of the theory it is possible to write all possible second order transport 
coefficients appear in the expression of J^. However, like energy momentum tensor, from the 
expression of retarded Green's function it is not possible to compute all the transport coefficients 
that appear in different order of derivative expansion. 

In this section we study the flow equation of retarded Green's function of boundary R- 
current. Unfortunately we find it difficult to solve the flow equation analytically to extract any 
transport coefficient. We present our calculation how to write the flow equation for retarded 
Green's function of i?-currents in presence of generic higher derivative terms in bulk Lagrangian 
and some numerical results. 

We start with Einstein-Maxwell action 

s = i6k/ dW ^(* +12 -H' (209) 



Solution is given by equation (196) with a' = 0. The temperature of the black hole is given by, 

A. 



T=fil-^\ (210) 



and the chemical potential is given by, 



H=—^. (211) 
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For technical advantage we write the metric and gauge field in a different coordinate. We 
change the radial coordinate r — >■ In this coordinate the metric and gauge field is given by, 

r 4r 2 U(r) r 

A t (r) = E(r) (212) 



where, 



2 r 2 \/3Q 
f/(r) = (l-r)(l + r- V). E(r) = ^(l-r). (213) 



We turn on small fluctuations for x component of gauge fields. Since the A t component 
of the bulk vector is non-vanishing in this background, the perturbations A x can couple to 
the tx component of graviton Therefore we also need to consider small metric fluctuations for 
components gt x . Writing them in momentum space 

Mr,x) = J ^i eik - X Mr,k) 

9t(r,x) = J ^ e *-*$(r,fc). (214) 

However, there exists a constraint relation between A x and g tx . We use this relation to replace 
gt x from equation of motion of A x . 

The on-shell action for gauge field fluctuations are given by, 



Sa = 



d 4 k 



167rG 5 J (2tt) 4 
The current corresponding to A x fluctuation is given by, 



(215) 



Ur,k) = jJ^ = - 2 r 2 U(r)A> x (r,k). (216) 



The equation of motion for A x (r, k) is given by, 

(U(r)K(r, k))' = _ t^j Ax (r, k) - E>(r)<f>\r, k) 



' ' J Vl,.(r./.-! • r/ ' ;,( ; !i .l,i r./.-). (217) 



4rQ \rU(r) r 

Using the constraint relation (coming from rx component of Einstein equations), 

^r,k) = - rE ' ir) i^ k \ (218) 
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and the equations of motion one gets, 

r( b\ - -( u2 _ t. 
Ar > } ~ 2 \rU{r) r 

Next we define a response function 



a(r, k) 



A x (r, k)-2rE'(r) 2 A x (r, k). 



J x (r, k) 



iuA x (r, k) 



(219) 



(220) 



Taking the derivative with respect to r and using the equation of motion we find the flow is 
given by, 



er'(r, k) 



2r U(r) _ 



a(r,kf 



r o ( q 2 U(r)\ 4r rE'(r) 2 U(r) 



(221) 



From the regularity of the response function at the horizon we find the boundary condition is 
given by, 



<T(l,k)' 



(222) 



With this boundary condition one can integrate this nonlinear equation to find finite fre- 
quency response of boundary Green's function. 




0.2 0.4 0.6 0.8 1.0 



Figure 3: Flow of R current correlation function. 
In presence of generic higher derivative terms in the Lagrangian the on-shell action for 



fluctuation A x may not have canonical form like (215). In that case one has to write an 



effective action for A x like transverse graviton. The effective will have the same form as (215) 
only the coefficients will depend on coupling constant of higher derivative terms. 



We conclude this section by presenting some numerical solutions of the flow equation (221) 
for both non-extremal extremal black holes in fig. [3] 



11 Conclusion and future directions 

In this review we have discussed the hydrodynamic behavior of strongly coupled systems in 
the context of the gauge/gravity correspondence. We have learnt how we can apply the holo- 
graphic methods to compute different transport coefficients of boundary fluid systems. We have 
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also elaborated the effect of higher derivative terms, which come from string theory, on these 
coefficients. 

First, we have shown how the shear- viscosity coefficient can be viewed as the effective 
coupling of transverse graviton in a two derivative theory. Next, we generalized this statement 
to any arbitrary higher derivative theory. We have developed a procedure to construct an 
effective action for transverse graviton in canonical form in presence of any higher derivative 
terms in bulk and showed that the horizon value of the effective coupling obtained from the 
effective action gives the shear viscosity coefficient of boundary fluid. Our results are valid up to 
first order in \i (coefficient of higher derivative term). We discussed two non trivial examples to 
check the method. We have considered four derivative and eight derivative (Weyl 4 ) Lagrangian 
and calculated the correction to the shear viscosity using our method. We found complete 
agreement between our result and the results obtained using other methods. 

In section [6] we have discussed that the shear viscosity to entropy density ratio is controlled 
by the near horizon geometry of the dual black hole spacetime. Therefore the knowledge of 
the near-horizon behavior of different bulk fields are enough to compute the ratio. That is why 
the computations are not very complicated even in the presence of higher derivative terms. 
We would like to emphasize that there is no ambiguity in defining the overall coefficient of the 
effective action and the results are consistent. We apply this idea to compute rj/s in section 



6.3 and see that the ratio is controlled by the horizon values of the scalars and so, in the non- 
extremal case, an operator deformation in the QFT will produce an interpolating non-trivial 
flow in which the moduli approach the (IR) black hole horizon. In the extremal case, though, 
the horizon moduli values are fixed and so the shear viscosity to entropy density ratio does 
not depend on the asymptotic values of the scalars. Therefore, QFTs with different UV fixed 
points can flow to the same IR fixed point. 

One important point to mention here is that the radial independence of retarded Green's 
function (response function) depends on the massless properties of transverse graviton h xy 
at k — > limit. We have seen that the mass term is proportional to k 2 and it vanishes 
in low frequency limit. However, in general this observation is not true. If we have other 



matter fields in the Lagrangian (like magnetic field considered in section 6.1) then all the 
components of transverse graviton may not be massless. Please look at appendix[Ej The isotropy 
(SO(3) symmetry) in spatial direction is lost in presence of a magnetic field. Therefore though 
xy component of transverse graviton (h xy ) is massless but other components h xz , h yz are not 
decoupled form gauge field perturbations and hence not massless even in low frequency limit. 
Therefore the retarded Green's function (response function) corresponding to these massive 
fluctuations are not independent of radial direction. Hence, horizon value and asymptotic 
values are different in these case. As a result the isotropy of viscosity tensor is also lostp*} 

We have also studied the flow equations for retarded Green's function of boundary theory 
analytically and found higher order transport coefficients of the boundary plasma solving this 
equation. We have generalized the analysis for generic higher derivative gravity theory. The 
flow equation for Green's function is a first order non-linear differential equation of Riccati 
type. Because of its non-linear nature it is hard to solve this equation exactly. After a change 



2!) 



Since, [T xy , T xy ] correlator is different than [T XZ ,T XZ ] or [T yz ,T yz ] correlators. 



51 



of variable one can reduce this non-linear equation to a second order linear homogeneous dif- 
ferential equation. But to solve this we need to specify two boundary conditions. In this article 
we have dealt with the non-linear equation and specify the boundary condition at the horizon. 
Therefore the hydrodynamic characteristic of the field theory at UV fix point is determined by 
IR boundary condition. In this way of computing the transport coefficients has an advantage 
over usual Kubo approach. In Kubo approach, one has to first find the transverse graviton 
by solving a second order differential equation and then compute regarded Green's function. 
Instead, the flow equation is a first order differential equation (although non-linear). As we 
want a perturbative expansion of Green's function in powers of u and q the equation turns out 
to be a linear first order differential equation. Thus, technically, it is simpler to get results for 
causal hydrodynamics, particularly when the dual bulk theory is complicated. 

It would be interesting to use renormalization — group flow techniques to understand how 
the response function evolves as a function of the radial direction. One can also apply Blackfold 
techniques of [52] and the idea of Wilsonian RG flow explained in [53J. Blackfolds are higher 
dimensional black holes with two separate length scales (mass and angular-momentum) and 
their equations of motions are obtained from fluid correspondence (in leading Einstein gravity). 
Generalizing this technique in reverse, one should obtain equations of boundary plasma from the 
Blackfold equations in a generic gravity theory. It would be really nice to check this explicitly. 

There are other methods to compute first and second order transport coefficients holograph- 
ically. We have discussed two of them in little details in appendix [Gj In G.l we have briefed the 
method developed in [28] . This is an elegant method to understand the properties of boundary 
fluid from boosted black brane geometry. This method shows that, all the constituent equa- 
tions of fluid system are captured in the Einstein equations. In [39], this method has been 
extended to higher derivative gravity (Gauss-Bonnet gravity) and shown that the above men- 
tioned statement is true even in presence of higher derivative terms. Although this method 
captures almost all transport coefficients (except k, the second order transport coefficient) up 
to second order in derivative expansion, but it is technically very hard to extend this idea to 
study higher derivative corrections. One problem is that the exact form of boundary stress 
tensor is not known for any generic higher derivative gravity, for example Weyl 4 term. 

In G.2| we have discussed another method to compute first order transport coefficients 
developed in [M] . In [22] , the authors proposed a Wald like formula for shear viscosity coefficient 
in generic higher derivative gravity. Though the formula works for generic four derivative term 
in bulk action, but unfortunately it fails to produce the correct result for Wely 4 term [10J. 
Motivated by this formula the author in [51] proposed a new Wald like formula for the first 
order transport coefficients. The author has shown that the transport coefficients can simply 
be obtained by evaluating the residue of the pole of the quadratic action near the horizon. It 
only assumes the massless behavior of low energy perturbation and regularity of the field at 
horizon. The result is independent of the boundary terms and thus one can neglect them. 

We conclude this review with some discussion on the currents trends on this subject. 

In this article we have mainly considered the characteristic of conformal fluid which has 
zero bulk viscosity at first order. Bulk viscosity measures the reaction of fluid against any bulk 
stress. The holographic computation of bulk viscosity is also an interesting topic in this area 
of research. In [55] the author has proposed a holographic bound on bulk viscosity to entropy 
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density ratio of strongly coupled gauge theory plasma, 



2 



>=---<£, (223) 
s 2n \p J 

where c s is the speed of sound and p is the spatial dimension of the system. This bound is 
dynamical, i.e. the right hand side depends on temperature unlike the bound on shear viscosity 
to entropy density ratio. They observed that the bound is saturated by the p + 1 space-time 
dimensional gauge theory plasma holographically dual to a stack of near-extremal flat Dp-branes 



at leading order (i.e. infinitely large 't Hooft coupling). The bulk viscosity bound 223 is also 
saturated in toroidal compactifications of conformal theories [55, 56J. It has been found that in 
various examples of string theory embedding of the holographic gauge-gravity correspondence 
the bound is saturated. For example, Af = 2* gauge theory plasmas, gauge theories with adjoint 
R charges, cascading gauge theories etc. There are further evidences in favor of saturation of 
this bound in [571 EE] for phenomenological models of fluid-gravity correspondence. However 
[59, 60] considered some other phenomenological models where the bound is violated. Recently 
in [SI] the author addressed the question if the violation of bulk viscosity bound is only limited 
to the phenomenological models of gauge/gravity correspondence. They considered strongly 
coupled M = 4 supersymmetric Yang-Mills plasma compactified on a two-manifold of constant 
curvature and found that the resulting (1 + l)-dimensional hydrodynamic system can have 
bulk viscosity coefficients which violates the bound if the curvature of the compact manifold is 
negative. It would be also interesting to check this bound explicitly in presence of some higher 
derivative terms in bulk. 

Non-relativistic generalization of the original AdS / CFT correspondence has opened up new 
directions in current research after a series of beautiful experiments on cold atoms at unitarity. 
The correspondence has been extended to explore the holographic duals of strongly coupled 
non-relativistic conformal field theories. Cold atoms are system of fermions interacting through 
a short-range potential which can be fine-tuned to obtain a massless bound state. Since the 
theory is scale invariant, it can be described as a non-relativistic conformal field theory with 
Schrodinger symmetry. This symmetry consists of the usual Galilean invariance, the scaling 
symmetry as well as the particle number symmetry. There has been some activity along this 
direction where solution generating techniques have been used to obtain bulk geometries with 
asymptotic Schrodinger symmetry. 

At the first order in derivative expansion, we have one interesting transport coefficient, which 
is shear viscosity coefficient. It turns out that for non-relativistic fluid the shear viscosity to 
entropy density ratio saturates the KSS bound [62J. However, [62] has also computed other 
transport coefficients like, thermal conductivity kinematic viscosity, thermal diffusivity, Prandtl 
number holographically, using the dual geometry proposed in [H31 EH ES] • 

As mentioned earlier, fluid dynamics in general can be thought of as an effective field 
theory with an infinite number of irrelevant terms obtained as usual in a derivative expansion. 
In relativistic case consideration of second order hydrodynamics is essential because the first 
order formalism is inconsistent with causality issues. However in non-relativistic setup one is 
not forced to consider the second order (beyond ideal and viscous terms) terms in stress tensor. 
But it is interesting to study those terms if we want to view hydrodynamics as a derivative 
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expansion of its local variables. 

One interesting problem in this direction is to study the light-cone reduction of relativistic 
conformal/non-conformal hydrodynamic stress tensor [28j |29| l30~t |4*6| 1^5] . In second order we 
encounter different transport coefficients in case of relativistic theory (see section [7]). It is 
possible to reduce the theory along one of the light-cone coordinates and understand how the 
higher dimensional transport coefficients descend down to new transport coefficients for the 
lower dimensional non-relativistic theory. The reduction of first order stress tensor has been 
done in [63J. Following the same procedure one can write the second order corrections to the 
stress tensor of a charged fluid with Schrodinger symmetry. Once we have all possible second 
order transport coefficient of a non-relativistic fluid we can evaluate their holographic values 
using the non-relativistic gauge/gravity correspondence. 

To compute the second order non-relativistic charged fluid's transport coefficients holo- 
graphically, one has to take the bulk spacetime to be five dimensional Reissner-Nordstrom 
black hole (embedded in 10 dimensional geometry) and perform TsT transformation to get 
asymptotically Schrodinger spacetime. Then one should turn on hydrodynamic fluctuations 
in this background and study the flow of corresponding response function (retarded Green's 
function) like |31j . Evaluating the response function at asymptotic boundary it is possible to 
compute the non-relativistic second order transport coefficients of 2 (spatial)-dimensional dual 
fluid system. However, in this way of computing the transport coefficients (Kubo method) 
may not capture all possible coefficient (like relativistic case). Applying the method described 
in [28] ([63] has already applied this method to evaluate first order transport coefficient for 
non-relativistic fluid) would be a good idea to compute the second order transport coefficients 
holographically. 
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A Second order transport coefficients from Kubo for- 
mula 



In this appendix we re-derive second order transport coefficients k and t w from usual Kubo 
approach. Let us now consider the action with solution given in section [7j 

' 1 - ^[R+12]. (224) 



S 



167tGh 



(TXy 



For well defined variation of this action we need to add a Gibbons-Hawking boundary term. 
Also, requiring the on-shell action being finite at boundary, we have to add counter-terms 
following usual approach of holographic renormalization. They are as follows: 

1 



Sgh 

ScT 



8ttGs 



d x\J— 7/C, 



16nG 



(225) 



(226) 



where 7 and 1Z are boundary metric and Ricci scalar (constructed out of 7) respectively. 
We consider the following metric perturbation, 

Qxy = 9xy + h X y (j", x) 

where e is an order counting parameter. We are interested in quadratic on-shell action for this 
transverse graviton $(r, x). Let us define the Fourier transform, 
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' ' - ik -*^(- ^ (227) 



<f>(r, k) 



(2tt)< 



$(r, x) 



and k = {—u, k}. Substituting this fluctuation in action (224), we get (236) 



HO 



Now, we rewrite 



this action (224) as equation of motion piece (which vanishes on-shell) and boundary term 
Thus, on-shell S, Sgh, Sct become, 

d 4 k 



31 



s 

Sgh 
Sct 



1 



16ttG 5 
1 

8ttG 5 
1 

167rG fi 



=6 (2tt)< 

d A k 

d 4 k 



-M<P(r,k)), 



^gh 



(<f>(r,k)), 



(2k)' 



L CT (<f>(r,k)), 



(228) 



where we define, 



L GH {<j>(r,k)) 
L CT (<f>{r,k)) 



a2(r)<f) , (r, k)<p*{r, k) + a4 ( r ) ^'^ ^ fc )0*( r , k), (229) 



2(^l)0(r, k)F(r, k) + 2(jg2)(f> , (r, k)<f>*(r, k), 
(cq + cxu? + c 2 q 2 )4>(r, k)(f)*(r, k). 



(230) 
(231) 



30 ' 0,3,0,6 are zero here 

31 We ignore contribution from horizon as |50j 
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The coefficients 02,04,05 are given in (237) and other coefficients are, 



2 - r 2 1 - r 2 Vl-r 2 1 y/1 - r 2 . 

9i = — 2~, 92 = -2 , c = -3 — , ci = - — c 2 = — (232 

r 4ry 1 — r 2 4r 

The retarded Green's function is defined at asymptotic infinity as, 

G n = 2 lim (L + L GH + LcT) \ on -s h eU _ 



Here, (f>o(k) is the boundary value of the fluctuation ( |227 ). The retarded Green's function 



G R is a function of boundary momenta k^ = (u,0,0,q). Now the leading action and the 
Gibbons-Hawking action get divergences from (f>'(r, k)<p*(r, k) and <p(r)<p*(r,k) parts. Both 
these divergences get canceled by the counter-term action which is always proportional to only 
4>(r)(f)*(r,k). In this case of leading Einstein's gravity, it is even more simplified. Divergences 
coming from <ft(r)<f)*(r, k) piece of leading and Gibbons-Hawking action gets canceled by mo- 
mentum independent piece of Counter-term action. It turns out that there is a cancellation 
among the corresponding coefficients as, 

lim(i(o4(r) - a5'(r)) + 2gl + c ) = \, (234) 
1 — >o z z 

i.e. the final contribution from <j>(r)<j)*(r,k) piece is only a finite number |. As the graviton 
fluctuation <f>(r,u,k) = 4>o(l + F(r,co,k)) and moreover \im r ^ F(r,co, k) = 0, we see that the 
cf){r)(j)*{r) term above only contribute to pressure (the oj independent piece of G R ). It would 
never contribute to any transport coefficient. 

Also, the divergences coming from <ft'(r, k)4>*(r, k) piece of original and Gibbons- Hawking 
action, get canceled with the piece of the counter-term proportional to uj 2 ,q 2 (ci and c 2 are 
purely divergent at boundary). Here, the situation is more subtle, as there is no cancella- 
tion among the coefficients. One actually needs to put the solution of <p(r,u,k) to see the 
cancellation. 

The overall lesson from this detailed analysis is that counter-term only cancel the UV 
divergences in usual holographic renormalization process and at most contribute to pressure of 
the boundary plasma. It has no effects on any transport coefficients. In [29] , the author have 
computed second order transport coefficients for the plasmas dual to leading Einstein's gravity 
following this usual approach. The results are as follows, 

TV = ^^, k = ± (235) 



These results match with the one we obtained in (125) by solving the flow equations 



B Equivalence of Boundary Terms 

In this appendix we will show explicitly why the transport coefficients computed from the 
original action and the effective action are same, even for any higher derivative theory. It was 
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already noticed [10], that the two would give same first-order transport coefficient rj with a 
suitable choice of the overall normalization constant. Here, we show that, not just the first 
order transport coefficients, rather any higher order transport coefficients computed from the 
original action and the effective action are same. 

We consider a general class of action for which appears when the higher derivative terms 
are made of different contraction of Ricci tensor, Riemann tensor, Weyl tensor, Ricci scalar 
etc. or their different powers. Since, all these tensors involve two derivatives of metric they can 
only have terms like d a d},§(r, x) and its lower derivatives. Therefor the most generic quadratic 
(in $(r, x), in linear response theory) action for this kind of higher derivative gravity has the 
following form (in momentum space|^ 

7-1 



<7 - -JL_ f d k ' 



al(r)(j)(r) + a2{r)(f){r) + a4(r)0(r)0'(r) + a a6(r)</>" (r)<f>' (r) 



+a' a3(r)0"(r) 2 + a5(r)0(r)0"(r) (236) 



where, 



al(r) = ^jf^ 8 + ol /2(r), a2(r) = -3r + \ + a' h2(r) 

6 4 
a 4(r) = — - - 2 + a' g2(r) a5(r) = -4r + - + a' j2(r) (237) 

and a3(r), a6(r),j2(r),g2(r), h2(r) and f2(r) depends on higher derivative terms in the action 
and hence are computed purely from the background solution with a' — > 0. Among these 
coefficients a3 is special, as, it couples to (J)" 2 . All four derivatives act on the graviton fluctuation 
and thus a3 only depends on background metric functions and there r-derivatives. It is easy to 
convince ourselves that a3 oc r(r 2 — l) 2 /(r, a'), where f(r, a') is a function that depend on the 
higher derivative terms and finite (constant or 0) at the boundary r — > 0. Now let us write the 
effective Lagrangian as follows, 



1 + aT f d 4 k 
16ttG 5 J W) 



S e ff = : - „ / T^rdr 



4r (r 2 - I) 2 0'(r) 2 - U W + ^ , + 



4r 2 (r 2 -1) 



(238) 



Demanding that the equation of motion (up to order a') of <ft derived from the original action 
and the above action are same we get, 



bl( r ) = _^[(-4r 3 - 12r + w 2 )a3(r) + (r 2 - l)(2nr A - a6'(r)r 3 - 4/tr 2 + 2a3'(r)r 2 



1 

2r (r 2 - 1) 

+2(r 2 - l)h2{r)r - 2(r 2 - l)j'2(r)r + a&(r)r + 2k + 2a3'{r))} (239) 



:',2 



In all the expressions we have omitted k dependence of cf). 
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62 (r) 



u A + 144rV)a3(r) + 4(r 2 - 1) ( - 4r 2 /2(r)(r 2 - 1) 



16r 2 ( r 2 _ X )4 

+ ((w 2 k - 2r 2 (r 2 - l)j2" (r))(r 2 - 1) + 2r 2 #2'(r)(r 2 - l) 2 + rco 2 a3" (r))(r 2 - 1) 
+ (1 - llr 2 )w 2 a3 / (r 



(240) 



The boundary terms coming from the original action (after adding Gibbons- Hawking boundary 
terms) are given by 



16ttG p 



d A k 



+ 0(r) 2 + rcf) (r)(j)(r) — 



(f>'(r)(f)(r) 



/1 1 nf,'(r) 

+ a '^-^2(r)0(r) 2 - -j2'(r)0(r) 2 + (W(r) - j2(r) - -Aiy( r )0( r ) 

a3'(r) (0(r)w 2 + 4 (r 4 - 1) <//(r)) 0(r) a3(r)(6r0(r)0'(^ 2 ) 



4r( r 2_i) 2 4r(r 2 -l) 3 
a3(r) ((r 2 - 1) (8r 3 + 24r - to 2 ) <p'{r)) 4>{r) a3{r)(j)'{r) (0(r)w 2 + 4 (r 4 - 1) 0'(r)) 



— a3(r)<//(r) 



4r (r 2 - l) 3 
0(r)u; 2 (r 4 -l)0'(r) 



4r (r 2 - l) 2 



2r (r 2 — l) 2 r (r 2 — l) 2 



(241) 



And the boundary terms coming from the effective action are given by, 



1 



d A k 



a 



r- - 0(r)0'(r) + - 2 
r) 2r(r 2 -l) 



-)(2r (r 2 -l) 3 + (-a6'(r)r 3 



167TG5 J (2?r) 

+2a3'(r)r 2 + 2 (r 2 - l) /i2(r)r - 2 (r 2 - l) j2(r)r + aQ'(r)r + 2a3'(r)) (r 2 - l) 



+ (-4r 3 - 12r + co 2 ) a3(r))0'(r)J 



(242) 



Now, it is interesting to compute the difference between these two boundary terms and the 
result if^ 1 



qB qB 



16ttG f 



d A k 

a3'{r)u 2 <f){r) 2 a3(r) (6rw 2 ) <p{rf 



l '' r + 0(r) 2 + a'( ^2(r)0(r) 2 - ij2'(r)0(r) 2 



4r (r 2 - 1) 



4r (r 2 - 1)' 



(243) 



The term proportional to a3 in the parenthesis of (243) vanishes at the boundary whereas the 



term proportional to a3' gives a pure UV divergent piece and a vanishing piece, due to the 
property of a3 mentioned above. This is true irrespective of the choice of the higher derivative 
terms. Thus, we see that the two boundary terms differ only by terms which are either purely 



33 There was a sign error in [10] 

34 It has been shown in [TU] that T = 0. 
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divergent or of the form g(r)<j) 2 , where g(r) is any function of r. The divergent terms would get 
canceled once appropriate boundary terms are added (which has been discussed in sections Pf 



and 7.3). The terms proportional to <p can only contribute to pressure of the boundary theory 



and are not important for the computation of transport coefficients of the boundary plasma. 
Thus we see that, it is obvious that the transport coefficients coming from the original action 
and the boundary action are same. 

Here we have considered only gravity theory. A more rigorous proof is required for 
theories involving covariant derivatives of curvature tensors and scalars. 

C Functions appeared in four derivative action 



A$ B (r,k) 



q 2 (r 2 — 1) + u 2 n 



+ 



g 2 (2ft(13r 2 - 3rco 2 - 13) + 130(r 2 - l)ft 



4 r 2 ( r 2 _ 12r 2 ( r 2 _ i 

+ (-36r 4 + 25r 2 + ll)ft) + oo 2 ((6r 2 - ll)ft + 130ft + 26ft) 



and 



A^ B (r,k) = r 



1 (r 2 - 1) ((18r 2 - 13) ft + 110ft + 22ft) a' 



3r 



(244) 



(245) 



D Coefficients in expansion (84) at 0(r — ru) 2 

For completeness, in this appendix we present the other coefficients that appear in the near 



„,„ , 5B 2 , 7q 2 \ , 3^ 2 



horizon expansion of different fields in eq. (|84|). They are 

u 2 = e 

1 / 

v 2 



3r 4 + 3 )^4^' h 



ziri - 2 



L_ L-^h ^245 4 C 2 e a * h + 24 (^>B 2 + 2g 2 ^ e 3a<t,h - 72B 2 ( 2 
(qB 4 [a 2 -A- B 2 q 2 (da 2 + 8j + q 4 ^3a 2 - a\\ e 4a ^ - U4e 2a ^ \ 



288vr 
+ 



-4acj> h 



W 2 



288tt 2 T 2 



36fi 4 C 4 - 24S 2 C 2 [B 2 + q 2 ) e 3a0h + 48 ( q 2 - 2B 2 ) e Wh + 288B 2 ( 2 e 2a ' t ' h 



3B 4 (a 2 - 4J - 8B 2 q 2 + q 4 (^a 2 -4^ e 6a<f>h - lUe 4a<f "^j 



22 



B( (B 2 (9C 2 e- 3Q ^ + 1) + 5 (q 2 - 6e" Q ^)) 
67rT 



Bq(e~ 3a ^ ((3a 2 + 4) (B 2 - q 2 ) e 3a * h - 12B 2 ( 2 + 24e 2a< M 
Pi = 96^ • (246) 
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E Decoupling of h xy mode 

In what follows, we provide a detailed derivation of the decoupling of the dual gravitational 
mode (4.1). We have explicitly checked that the h xy = e ltu)+2V ^e&(r) mode does not couple 
with any other field when the momentum vanishes. 



For two derivative gravity theory, this can be easily seen from the equations of motion (66 



67, 68). However, we are interested in the most general four-derivative action (91). In this 
case, instead of writing the equations of motion in the presence of higher derivative terms, we 
will explicitly compute the action up to order e 2 . In this way, it can be explicitly checked that 
there is no coupling between h xy and the other fields. 

Let us turn on the following perturbations of the metric 

9 a/3 = 9^1 + e haf3 



( 



V 

and gauge gauge fields 

Atu) , 





e ituJ e& 


e itw e£ 3 


e itw e£ 4 


e ituJ e& 




e 2W Z 2 - U 


e itw eT 





e 2W Z + e itu 




e <tw eT 


e 2V 


e itu+2V e Q 


e itu ex 







e itoj+2V e Q 


e 2V 





e itul e& 


e 2W Z + e ituJ ev 


e itu ex 





e 2W 



ev 



(247) 



By 



Bx 



e ltw ea r (r), e ltU} ea t (r) - E(r), — + e ltU} ea x (r), + e lt "ea y (r),e ttuJ ea z (r) + P(r) 



(248) 



Here g^l and are the background metric (|71|) and the background gauge field; e dependent 
terms are the perturbations. 



Using these field excitations, we can now compute the action _ The result is complicated, 
but for our purpose it is enough to pick up the $(r)-dependent terms — for concreteness, let 
us write the $(r)-dependent part of the action: 

35 We use the Mathematics notebook for this computation. We emphasize that our system is symmetric in 
x— and y— directions. 
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s. 



4C/(r 



-((e _2y(r)_M/(r) (-f/(r)e 2W(r) (4i? 2 e a</,(r) + AQ(r) 2 e 4V{r)+aip ^ r) 



-2e 4V ^U"(r) - 4e 4V »f/ / (r) (2V'(r) + W'(r)) + e 4y W+ 2W Wz'(r) 2 + 24e 4V ») 
+2U(r) 2 e 4V{r) (2P'(r) 2 e a ^ r) + e 2W{r) (AV"(r) + AV\r)W\r) + QV'(r) 2 

+2W"(r) + 2W\r) 2 - ^'(r) 2 )) + Uu 2 e 4V ^ +2W ^))) 
+$(r)$'(r) ( 2e 2y M +H/ M ([/'(r) + U(r) (3V'(r) + W"(r))) 

-(2a'e H/ M- 2y M(-2f/(r)\/ / (r)(f/(r) (B 2 c 2 + Cl e 4V ^ +2W ^ Z'(r) 2 ) - 2 Cl f/(r) 2 e 4V ^ 
(2V"(r) + 3\/'(r) 2 + W^'(r) 2 ) + 2 Cl u 2 e 4V{r) ) + 2c 1 f/(r)e 4V ' (r) f/ / (?f^ / M 

+c 1 e 4y Mf/ , (r) (2U(r) 2 (V"(r) + 3V'(r) 2 ) + uj 2 )))/U(r)) 



+$(r)$"(r 
+$'(r) 2 



(2f/(r)e 2V > )+H/(r) - A Cl a'U{r)e 2V ^ +w ^{U'{r)V{r) + 2C/(r) (V"(r) + V'{r) 2 ))) 



a'e w ^- 2V ^(U(r)(B 2 c 2 + 2c 1 e 4y ('V(r)\/ / (r) - Cl e 4V ^ +2W ^Z'(r) 



+ Cl e 4V ^ (U'(r) 2 + Au 2 ) + 2 Cl U(r) 2 e 4V ^ (-2V"(r) + V\rf + W\r) 2 )) + ^[/( r )e 2y ^ + ^) 



+$'(r)$"(r) 



(2c 1 a / f/(r)e 2V 'M +H/ M (C/'(r) + AU{r)V'{r))) 



+ $"(r) 2 



2c 1 a / f/(r) 2 e 2y W +w W 



(249) 



Since there is no mixing between <3> and the other modes in the zero momentum limit, this 
mode remains massless and a computation of the related Green's function from the near-horizon 
data is still possible. 

However, the other metric fluctuations namely h xz or h yz are not decoupled from the gauge 
field perturbations in the zero momentum limit (even in the absence of higher derivative terms) 
- these coupled terms are proportional to B or ( (CS term). Therefore, these modes are not 
massless in the zero momentum limit and the near horizon geometry is not sufficient to compute 
the related two point correlation functions, i.e. < [T XZ ,T XZ ] > or < [T yz ,T yz ] > . 
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F Functions appeared in Higher-derivative Charged black- 
hole action 



CB 



and 
A C X B 



r 5 (u 2 - q 2 ) + q\%r da' (llr 8 (to 2 - q 2 ) - r 4 r 4 (25g 2 + Qu 2 ) + 36<fr«) 



+ 



2(r 4 - r 4 ) 

Q 2 



6r 3 (r 4 



'0/ 



a' 



rV + — ( d[28q 2 r®r 2 - r 8 (36g 2 + 127^ 



2r 2 (r 2 — r 2 )(r 2 + r 2 ) 2 
-4r 2 r 6 (7g 2 - 3w 2 ) - 8g 2 r 8 + 4r 4 r 4 (llg 2 + 3w 2 )] - 24c 2 [r 2 r 6 (2g 2 - 3w 2 ) 



+ r 4 r 4 (2g 2 - 3tu 2 ) - 2q 2 r* + 4r V 1 



0(Q 4 



2q 2 r%r 2 
(250) 



1/4 5^ ^a' (13r 4 - 18r 4 ) (r 4 - r 4 ) Q 2 (r 2 -r 2 ) 
2 l rr o r ) —» + 



6r 3 



2rr 2 



a 



+c 1 (101r 4 - 156r 2 r 2 - UOr, 



+ 0(Q 4 ). 



1 - — (24c 2 (4r 4 -3r 2 r 2 -3r 4 ) 



(251) 



G Other methods of holographic hydrodynamics 

In this appendix we outline two other methods of computing various transport coefficients very 
briefly. First we discuss about the boosted black brane method by [28J. Next, we talk about 
the pole methods by [54"] . 

G.l Boosted black branes and fluid dynamics 

In this section we will briefly sketch the working procedure of [2B]- For detailed discussion 
readers are referred to the original paper. We will also skip the technical details in this section. 



Consider the Einstein-Hilbert action with negative cosmological constant 



16nG 5 

where L is the radius of AdS space. 
The equation of motions are given by' 31 ' 



d^x^f^g R + 



12 



(252) 



Emn — Rmn — ^Rg ain - 



2 9mn — 0. 



(253) 



36 _Af 



x = {v, r, x}. 



62 



There exists a class of solutions to these equations of motion given by the "boosted black 
branes' : 



37 



^.2 ^,2 

ds 2 = —lundx^dr — —f{br)u^u u dx^dx v + — {u^u v + r]^) dx^dx v (254) 

with, 

and Uj = 7/3j (255) 



where, 7 = 1/ y 1 — f3 2 . 

Putting the values of u^s the metric can also be written as, 

^.2 ^,2 

ds 2 = 2jdvdr - —j 2 f(br)dv 2 + —dx l dx l 

L 2 L 2 
2 2 

+ T^(7 2 - l ) dv2 ~ ^l^dx l dr - 2 7 — 1 2 (l - f(br))f3 i dx i dv 
L 2 L 2 
2 

+ ^7 2 (1-/(^))^W. (256) 
The solution is parametrized by four constant parameters b and /3,'s. 

The black brane horizon is located at r# = 1/b and the temperature of this black brane 
is given by, 

T = (257 > 



Consider the metric (256 ) and replace the constant parameters b and /3j's by slowly varying 



functions 6(x M ) and /3j(x M )'s of boundary coordinates 

^2 ^,2 

(is 2 = 2jdvdr - —j 2 f(b(x a )r)dv 2 + —dx i dx i 

L 2 L 2 
2 2 

+ ^-( 7 2 - l)dv 2 - 2 1 p i {x a )dx l dr - 2^(1 - f^x^^x^dx'dv 
+ ^(l-fibix^rWiix^frix^dx'dxJ. (258) 
We will call this metric (b(x a ) , j3i(x a )) . 



In general the metric (258) is not a solution to Einstein equations unless one adds some 
corrections to the metric and also the parameters b(x a ), (3i(x a ) satisfy some set of equa- 
tions, which turn out to be the equations of boundary fluid mechanics. 



37 m 



x^ = {v, x}. 
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• Write the parameters b(x a ) and (3i(x a ) and the metric as a derivative expansion of the 
parameters. Up to first order in derivative expansion, 

g = g i0 \b(x a ), + eg m (b(x a ),Pi(x a )), (259) 

b( x a ) = b (0 \x a ) (260) 

and 

P l {x a ) = pf\x a ) (261) 

where e is a dimension less parameter whose power counts the number of (boundary) spacetime 
derivatives acting on the parameters. Since b^ l \x a ) and f3 < f'\x a ) do not enter in to the 
first order equation of motions, we have kept the expansion for b and /Vs up to leading 
order. 

• In general one can write the metric and parameters as power series of e. Then plug the 
metric in Einstein equations and solve the metric and the parameters order by order (in 
e). For example in our case since we are interested up to first order, we will plug the 
metric in Einstein equations and solve for g^ 1 ' and the constraint equations imply some 
relations between the zero" 4 order parameters. We will work in a particular gauge, 

Tr((0< o >)-y i} ) = 0. (262) 



After finding the metric with first order fluctuations one can find the boundary stress 
tensor (using the definition given in [HS])- The form of the boundary stress up to first 
order in derivative expansion is given by, 



16ttG 5 T^ = T 4 tt 4 L 3 (4tvi„ + rj^) - 2T 3 7r 3 L 3 a 



ilVl 



where is given by, 



and P^ u = u^u u + rj^ u . 



cr„„ = P^pPd^up) - ^P^ v d a u c 



(263) 
(264) 



One can follow this procedure to second order in derivative expansion. We are not discussing 
that lengthy algebra here. Rather we are going to present the final answer for stress tensor up 
to second order in derivative expansion. 



3 uv 



=(vr T) 4 (tT + 4 u»u v ) - 2 (tt T f a 

+ (ttT) 2 ( (In 2) TZ + 2 T% + (2 - In 2) 



I 2c I 1 2d 2e 



(265) 
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where 



= p^p^ d [a u p) - \ P^ d a u a 
T% = a^a u a - l -P^a^a aP 

Tg = d a u a (266) 
= Vu^ Vu v - i P^ Vu a Vu a 

Tg = P» a P^ V (d {a u p) ) - ± P a/S V (d aU(S ) 

Our conventions are eoi23 = — e 0123 = 1 and T> = u a d a and the brackets () around the indices 
to denote symmetrization, i.e. , a^ a b^' = {a a b^ + a l3 b a )/2. 

From this expression of second order stress tensor, one can read-off the transport coefficients, 

, T o^ 2 -ln2 , 2ri , 2ri In 2 , „„, 

r/ = -iV 2 T 3 , r u = — , A 1 = ^, A 2 = ^— , A 3 = 0. 267 

8 7X1 71 1 71 1 

As the coefficient k does not enter the equations of fluid dynamics in flat space (see equation 
(105)), this analysis leaves this coefficient undetermined. 



G.2 Pole Method 

The method is very well explained in the original paper [53]. Here we outline the working 
formula in brief. Let us consider a Lagrangian containing arbitrary higher-derivative terms. As 
usual, we treat the higher-derivative terms perturbatively small and the equations of motions 
can still be treated as quadratic. The most generic form of the quadratic action in Fourier 
space (at zero spatial momenta) is given as follows: 

S = J lifzldx 1 J ^(S z + S t + S B [<P(z,x\u)]), (268) 

Where S z is radial action, containing at least two radial (z) derivative and no time derivative. 
St contains time derivatives and are proportional to u 2 . Sb is the boundary Lagrangian. The 
only constraint that has been put to write the above Lagrangian density is that the field is 
massless. Studying the generic structure of the boundary terms, it can logically be shown that 
these boundary terms do not contribute to the computation of first order transport coefficients. 

The next step is to compute the quadratic action near the horizon. Now, by demanding the 
regularity at the horizon, it can be shown that the metric perturbation <fi has following solution 
(independent of the details of the action), 

(psoi = 0o exp -i^rp log z, (269) 
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where, T is the temperature of the background. The Lagrangian density computed on this 
solution will always have a pole and other regular terms in the radial variable z. The first order 
transport coefficient is then given by a simple formula as, 

„ ^Residuez = 0C-? , _ . 

S = 8vrT - 270 

ur 

where, £2 is the quadratic Lagrangian density evaluated at $ so \. This method is really simple 
and easily usable for theories with arbitrary higher- derivative couplings (powers of curvature 
tensors as well as their covariant derivatives) . But its extension to higher order transport 
coefficient is rather unclear. 
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